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1.1  General 

The  German  Experimental  Seismic  System  (GERESS)  has  been  operated 
by  Ruhr-University  Bochum.  The  GERESS  array  is  located  in  southeastern 
Germany  (Bavarian  Fmest)  near  the  border  to  Austria  and  the  Czech  Republic. 
The  array  consists  of  25  stations  with  vertical-component  short-period  Teledyne 
Geotech  GS-13  instruments  sampled  at  40  Hz.  Four  of  the  sites  include  hor¬ 
izontal  conqxMient  instruments.  At  the  key  statkm  of  the  array,  GEC!2 
(48.84511  N,  13.70156  E,  1132  m),  a  STS-2  is  sampled  at  80  Hz  (broad  band 
element).  The  geometry  of  the  array  is  based  on  concentric  rings  providing  an 
overall  tyrerture  of  about  4  km  (Harjes,  1990).  GERESS  was  installed  in  1990 
atul  the  number  of  recording  elements  increased  gradually  over  that  year.  The 
array  became  fully  operational  in  January  1991.  Data  from  the  array  have  been 
continuously  transmitted  to  the  NORSAR  data  center  in  Kjeller  and  to  the  Insti¬ 
tute  of  Geophysics,  Ruhr-University  Bochum  via  64-kbit  lines. 

This  report  summarizes  the  technical  perfonxiance  of  GERESS  and  the 
Bochum  data  center,  and  continues  similar  descriptions  given  previously  (lost, 
1991a:  1992a;  1993a,b).  This  summary  is  based  on  noonthly  status  repOTts  that 
are  available  upon  request 


1.2  GERESS  Data  Cento-  at  Ruhr-Univerrity  Bochum 

The  Institute  of  Ger^hysics  of  Ruhr-University  Bochum  has  operated  an 
experimental  on-line  processing  syston  for  GERESS  data  (Jost  1992b;  Harjes  et 
al.,  1993).  This  system  uses  software  develc^red  at  NORSAR  data  center  (Myk- 
keltveit  and  Bungum,  1984;  Fyen,  1987).  The  automatic  analysis  has  woiked 
very  stable  and  results  ate  available  for  all  GERESS  data  in  this  reporting 
period. 

The  on-line  processing  consists  of  3  steps:  detection,  fk-atudysis,  and  loca¬ 
tion.  The  first  stage  of  the  on-line  processing  accesses  data  in  30  second  seg¬ 
ments  and  runs  a  STA/LTA  detector.  The  detector  presently  recognizes  an 
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KEARSARGE  above,  causing  less  coupling  for  the  latter.  Gupta  et  al.  (1989) 
found  that  teleseismic  mi,  values  for  shots  below  the  water  table  are  0.3  •  0.5 
magnitude  units  larger  than  fOT  shots  above  the  water  table  at  a  fixed  yield. 
Following  these  authors,  we  assumed  that  the  water  table  influence  is  about  0.4 
magnitude  units,  scaling  our  observed  m),  =  5.4  (GERESS)  for  KEARSARGE 
to  5.8  for  a  regular  150  kt  NTS  shot  below  the  water  table.  We  obtained 

logic  (yield)  =  1-21  nib  (GERESS)  -  4.84  .  (1) 

The  resulting  yield  estimates  fiom  GERESS  are  given  in  Table  1-3.  For  com¬ 
parison,  magnitudes  and  yields  determined  at  GRF  by  Schlittenhardt  (1988b- 
1992)  are  included.  Using  Nuttli’s  (1986,  1987)  linear  relation  between  nob 
(Lg)  and  yield  and  taking  proper  care  of  the  magnitude  bias  between  P-wave 
and  Lg-wave  magnitudes,  we  obtained  furtiher  yield  estimates  firm  GERESS 
P-wave  magnitudes  (Table  2-3).  Recendy,  Goldstein  et  al.  (1994)  estimated  the 
yield  for  BRISTOL  as  10  kt  They  argued  that  their  value  could  be  30  %  lower 
due  to  differences  in  coupling.  We  estimated  for  BRISTOL  7  kt  from  (1),  and 
8  kt  using  Nuttli’s  relation.  Essentially,  the  yield  estimates  from  (1),  Nuttli’s 
relation,  and  GRF  are  not  significandy  different  considering  the  lack  of  detailed 
shot  information  (e.g.,  water  table,  ^all).  In  additicm,  effects  of  local  geology 
can  affect  teleseismic  P-wave  magnitudes  introducing  uncertainty  to  the  yield 
estimate  (McLaughlin  et  aL,  1987;  Cismier,  1987).  In  Table  2-3,  a  yield  esti¬ 
mation  was  calculated  to  all  events.  This  estimation  together  with  the 
automatic  detection  of  the  NPT  suggest  that  GERESS  can  detect  events  in  the 
order  of  1  kt  at  NTS. 

FurthermcHe,  Table  2-3  suggests  that  GERESS  can  observe  nuclear  tests  at 
NTS  clearly  below  a  yield  estimate  (tf  1  kt  Although  small  events  like  FLOY- 
DATA  or  GALENA  were  not  automatically  detected,  dieir  visible  waveforms 
can  be  used  to  confirm  a  location  hypothesis  within  the  fiamewotk  of  GSETT- 
3. 

In  order  to  identify  subgroups  of  the  Nevada  data  set  we  pettomed 
cross-correlations  of  the  t^timum  beams.  Hgure  2-7  shows  the  correlation 
matrix  obtained  by  using  the  maTinmm  cross-correlation  value  to  the  filter- 
band  0.5-2.5  Hz  and  die  event  numbers  firom  Table  2-3.  We  clearly  see  that 
the  small  events  (mb  <  4.3;  events  #  1  -  #  7)  do  not  have  a  correlation  larger 
than  0.54,  possibly  due  to  the  small  SNR.  On  the  other  hand,  the  cross¬ 
correlation  sequences  to  FLOYDATA,  AUSTIN,  DIVIDER,  and  GALENA 
indicate  that  the  observed  wavefonm  are  from  the  ccnresponding  events  at 
NTS.  For  the  larger  events  (mb  >  4.6;  events  #  8  -  #  19),  the  correlation 
reaches  values  of  up  to  0.91. 

We  further  performed  a  single  linkage  clustering  analysis  (Everitt,  1993). 
Hgure  2-8  shows  the  corresponding  dendrogram.  We  see  that  JUNCTION  (# 
17),  MONTELLO  (#  10),  HOUSTON  (#  14),  and  COMSTCXX  (#  13)  form  a 
cluster.  Also  KEARSARGE  (#  12)  and  BEXAR  (#  15)  tom  a  subgroup.  An 
explanation  to  the  latter  could  be  the  neaiiy  identical  shot-locations  (Fig.  2-9). 
AIm  HUNTERS  TROPHY  (#  6)  and  the  NPT  (#  5)  appttr  somewhat  related, 
although  the  significance  level  is  radter  low.  Obviously,  the  shots  at  Pahute 
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Table  2*5:  Magnitude  Yield  Relations  for  Tuamotu  Archipelago 


Yr 

Mo 

Day 

Name  /  Site 

“b 

GRF 

Yield 

GRF 

“b 

GERESS 

Yield 

DSIR 

No. 

1988 

05 

11 

Mururoa 

40 

5.3 

20 

~~r 

1988 

05 

25 

Mururoa 

65 

5.3 

83 

11 

1990 

06 

02 

Mururoa 

5.4 

40 

5.1 

30 

8 

1990 

06 

07 

Mururoa 

4.8 

10 

4.3 

3 

3 

1990 

06 

26 

Fangataufa 

5.7 

100 

5.5 

100 

12 

1990 

07 

04 

Mururoa 

5.2 

20 

4.7 

18 

5 

1990 

11 

14 

Fangataufa 

5.7 

120 

5.5 

117 

14  , 

1990 

11 

21 

Mururoa 

5.7 

35 

5.4 

36 

10 

1991 

05 

07 

Mururoa 

4.6 

4 

4.0 

1 

2 

1991 

05 

18 

Mururoa 

5.4 

30 

5.0 

16 

4 

1991 

05 

29 

Fangataufa 

5.8 

90 

5.4 

107 

13 

1991 

06 

14 

Mururoa 

5.5 

50 

5.1 

28 

7 

1991 

07 

05 

Mururoa 

<0.5 

<3.3 

<0.3 

1 

1991 

07 

15 

Mururoa 

5.6 

70 

5.2 

34 

9 

The  No.  column  indicates  a  sorting  of  the  events  according  to  yield  estimate  (DSIR). 
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Appendix  1-2:  GERESS  Instrument  Calibration 


Each  day,  the  GS-13’s  of  GERESS  are  calibrated  using  a  1  Hz  sinusoidal  input  voltage 
(GS-13  4500  mV)  fed  to  the  calibration  c(^. 


cal.  time  channel 


3K)0- 

3:01K» 

D6 

3:01 

3:01:30 

A2-m 

3.-02 

3.-02:» 

A2-ae,  B  -  ring 

3K)4 

3:04:50 

C1.C3 

2106 

306:30 

Dl-tz.  D2.  D3 

3:07 

307:30 

Dl-n,  D4^.  D7-ai 

3:08 

308:30 

Dl-ae,  D4-ie.  D7-sc,  C4,  C5,  C6.  C7 

3:10 

3:10:30 

D4-B.D3 

3:12 

3:12:30 

D7-n,  D8,  D9 

Calibration  on  May  26, 1994  (for  D6  from  May  8;  318273  cnts  are  expected) 


channel 

cnt 

deviaiion  in  % 

OEAl  iz 

319234 

0.308 

GEA2  az 

311215 

-2J12 

GEA3  sz 

316905 

-0/430 

GEBl~sz 

319442 

0.367 

GEB2~sz 

316260 

-0.632 

(XB3  fz 

380812 

1965 

GEB4  s 

320425 

0.676 

GEB3~sz 

308831 

-^96 

GECl  sz 

314933 

-1.05 

GEC2"sz 

312405 

-164 

GEC3  sz 

313323 

•0.926 

GEC4~sz 

311943 

-1.99 

GEC3  sz 

312094 

-1.94 

GEC6  sz 

315217 

-0.960 

GEC7  sz 

309903 

-263 

GEDl  sz 

297258 

660 

GED2  sz 

324986 

2.11 

GED3  sz 

318497 

0.070 

GED4  a 

300922 

-565 

(SD6_b 

317980 

-009 

GEOSja. 

312444 

-163 

GED7~8 

289007 

■920 

GED8~8 

316283 

-0625 

GED9~e 

320454 

0.683 

(XA2_a 

?»3890 

■0.749 

GEDl~n 

324707 

X02 

GEDl'a 

314142 

-160 

GED4~n 

312861 

-1.70 

OEl>4~a 

313429 

-162 

GEl>7~a 

309436 

-2.71 

GED7~a 

308842 

-2.96 
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derived  expressions  valid  to  second  order.  In  his  perturbative  treatment  the 
stiihiess  matrix  can  be  represented  as 

Cijkl  “  Cijtl  “I"  Cijkl  ^ijkl  (I^) 

for  long  wavelengths,  where  the  zero-order  term  pves  the  isotropic  reference 
model  specified  by  the  Lame  parameters  A  and  fi 

Cijkl°  =  ^Sij^kl  +  ft(SikSjl  +  SiiSjk)  (15) 

The  first  and  second-order  perturbation  terms  can  be  interpreted  as  first- 
order  effect  of  the  cracks  and  crack-crack  interaction.  For  the  crack  density 
e  =  Na^ /V  (N/V=number  of  cracks  per  volume,  a=mean  crack  radius)  they 
are 

A2  X(.  +  2fi) 

A2  A^  A(A  -1-  2/1) 

A(A-|-2/i)  A(A-(-2/i)  (A-|-2/x)2 

0. 


Cijkl  -  ^5 


X^q/iX  +  2fi)  X^q/(X  -I-  2fx)  Xq 

X^q/iX  +  2n)  X^q/(X  +  2ft)  Xq 

Xq  Xq  (A  -f  2fi)q 


with  q  =  15(  j)*  +  28(^)  -|-  28  and  x  =  Here  the  cracks  are 

perpendicular  to  xj.  D  is  the  diagonal  matrix  with  trace 


(17) 

aligned 


trace(D)  =  (un,Uii,«ii,U33,«33,0) 


Uim  depends  on  conditions  on  the  crack  face.  The  dependence  of  Cijki  on  the 
crack  parameters  is  only  through  D  and  is  further  given  by 


(19) 

and 

16  (,\-t-2p)  , 

3(3A  +  4^)"‘''"’' 

(20) 

with  k  and  m  given  as 

''  ird/i  X+ji) 

(21) 
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POLARIZATION  ATTRIBUTES  FOR 


square  amplitude  ot  the  three-component  seismogram  shows  the  mam 
phases  clearly.  Rectilinearity  and  component  ratios  allow  determina¬ 
tion  of  the  wave  type. 
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Figure  2-3:  Location  of  the  nuclear  test  sites  monitored  at  GERESS  (1988  - 
June  1994). 
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_  V  (A  +  2/i) 
irdfi  (3A  +  4/i) 

The  e-range  of  interest  in  this  study  is  between  0  and  0.1  .  For  such  low  val¬ 
ues  of  e,  the  second-order  approximation  is  considered  imcritical  (Crampin, 
1981).  d  is  the  (large)  aspect  ratio  of  the  cracks.  A'  and  fi'  are  the  Lame- 
parzoneters  within  the  cracks.  The  two  extreme  cases  with  particular  physical 
significance  are  1)  =  A'  =  0  and  2)  /i'  =  0  with  nonzero  A'. 

The  first  case  models  dry  cracks  (HCDl),  with  free  surfaces  at  the  crack 
faces.  In  this  case  k  and  m  vanish.  The  polarization  pattern  is  deter¬ 
mined  solely  by  (rj,  t)  =  (1.43,  2.13).  P-wave  phase  velocities  and  polar¬ 
izations  for  this  model  of  vertical  East- West  oriented  cracks  are  shown  in 
figure  3.9  in  lower  hemisphere  representations,  where  for  display  purposes 
the  symetry  axis  is  shown  pointing  upwards  in  the  paper  plane.  We  choose 
the  unconstrained  absolute  levels  of  Cn,  C44  and  cee  to  conform  to  model 
HCDl  (background  velocity  vp  =  5.8km/sec,  vs  =  3.349A:m/sec,  density 
p  =  2.6km/sec). 

This  choice  affects  only  the  P-velocity  pattern,  but  not  the  polarization 
pattern,  which  depends  only  on  rj  and  t.  Propagation  in  the  plane  of  aligned 
cracks  gives  high  velocity  and  the  P-wave  anisotropy  for  this  ideal  model  is 
24%.  The  right-hand  side  shows  polarization  anomalies,  namely  deviations 
of  the  polarization  from  propagation  directions.  In  general  the  polarization 
direction  deviates  into  the  direction  to  larger  phase  velocity,  diverging  at  the 
symmetry  axis  and  converging  to  the  symmetry  plane. 

For  the  second  case  (model  HCSl),  which  models  saturated  crzu^ks,  we 
obtain  figure  3.10.  The  parameters  are  the  same  as  in  figure  3.9,  but  with 
aspect  ratio  d=0.0001  and  A'=2.25  GPa.  The  polarization  pattern  is  then 
given  by  77=0.82  and  r=l.  Note  that  this  velocity  pattern  differs  substan- 
tizdly  from  the  dry-crack  model:  The  symmetry  axis  has  also  fast  velocity. 
This  gives  rise  to  t=1.  Clearly,  t  is  very  sensitive  in  distinguishing  dry 
crack  conditions  &om  saturated  ones.  The  difference  between  minimum  and 
maximum  velocity  with  3.5%  anisotropy  is  substantially  smaller  than  for  the 
dry-crack  model.  The  maximum  velocity  achieved  is  somewhat  higher  than 
for  the  dry-crack  model  though.  Polarization  anomalies  are  smaller  but  give 
complexer  pattern. 


3.5  MODEL  FITTING 

Observed  polarization  anomalies  can  be  related  to  the  parameters  tj  and 
T  describing  the  anisotropy.  The  inverse  problem  requires  solving  the  eigen¬ 
value  problem  (6).  In  this  paper,  this  is  done  by  nonlinear  parameter  fitting. 
As  shown  above,  parameters  in  the  inversion  are  rj  and  r,  but  also  two  angles 
0  and  4  describing  azimuth  and  dip  of  the  symmetry  axis.  In  general,  we 
wiU  not  be  able  to  determine  the  4  elastic  constants  independently.  Instead 


Observed  Polarization  Deviations 

(lower  hemisphere  projection:  0-50  deg) 
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Figure  3.3:  Lower  hemisphere  under  the  GERESS  axri^  (to  50”  inci¬ 
dence)  showing  event  polarization  anomalies  averaged  over  all  three- 
component  instruments.  Lower  hemisphere  locations  of  first  arriving 
phases  for  a  few  geographic  repons  are  shown.  Note  the  consistent  de¬ 
viation  to  Northern  directions  for  phases  steeply  incident  firom  North. 
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Rgure  1-2:  The  basic  menu  of  the  On-line  Processing  Display  Manager  at  the 
GERESS  data  center  Bochum 
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traces  start  10  seconds  before  the  theoretical  onset  time  calculated  from  the 
published  DOE  parameters  and  the  IASP91  travel  time  tables. 


02.34..45.100 


sorted  according  to  magnitude  determined  at  the  GERESS  site  (Table  2-9).  The 
trace  of  the  event  on  May  4,  1988  starts  10  seconds  before  the  measured  onset 
time,  the  other  traces  were  aligned  to  its  first  maximum  (the  start  time  of  each 
trace  is  indicated  to  the  right). 


66 


we  may  obtain  two  nonlinear  constraints  on  these  parameters  firom  rj  and 
r,  which  we  may  supplement  with  a  priori  information.  In  our  case,  such 
information  is  given  by  reasonable  assumptions  about  the  reference  velocity. 
In  the  following  we  will  use  reference  velocities  used  in  the  above  examples 
HCDl  and  HCSl.  This  enables  us  to  also  display  velocity  patterns  for  the 
inversion  results. 

For  each  event  we  have  an  observation  of  polarization  direction  (  ]  ^d 

propagation  direction  wish  to  compare  these  angular  differences 

^  with  predictions  from  models  depending  on  the  pa- 


f  Aei'{v,r,6,9)\ 
ting  the  misfit 

S(6,  r/,  r)  =  ^  ^2  cos~^  ^cos(A$i  —  A^,')cos(Ad,  -  A$/)^  j 


for  the  i-th  event. 


rameters  (  6,  #,  jj,  r  ),  which  we  call 
We  obtain  the  best  fitting  model  by  inspecting  the  misfit  on  the  sphere 

|2 


(23) 


which  is  the  analogon  of  the  standard  sum-of-squares  misfit  criterion  on  the 
sphere. 

The  models  HCDl  and  HCSl  have  shown  that  the  polarization  pattern 
may  vary  considerably,  depending  on  the  choice  of  q  and  t.  We  start  the 
discussion  of  inversion  results  by  studying  the  dependence  of  the  misfit  on 
symetry  axis  orientation  for  these  two  special  cases,  ij  and  t  are  fixed  to  (1.43, 
2.13)  in  the  first  case  of  dry  cracks  and  (0.82,  1.)  in  the  case  of  saturated 
cracks.  We  are  left  with  fitting  6  and  #.  The  nonlinear  search  may  use  a 
rotation  of  the  elastic  constant  matrix 


^ijkt  ■“  (2^) 

with  the  rotation  matrix  q.  While  this  allows  for  all  types  of  S3rmetry,  we 
may  more  efficiently  use  (9)  in  our  case  of  hexagonal  symetry.  Results  for  the 
dry-crack  case  are  given  in  figure  3.11  showing  the  achieved  misfit  S  given  in 
degree^  by  (23)  depending  on  the  orientation  of  the  symmetry  axis,  namely 
its  location  in  the  lower  hemisphere.  Axes  orientations  with  small  misfit  are 
near  the  drciimference.  The  smallest  achieved  misfit  is  135  degree^,  which  is 
substantially  larger  than  the  starting  misfit  of  68.8751  degree^.  The  simple 
dry-crack  model  is  clearly  incompatible  with  the  data.  On  the  other  hand, 
the  saturated-crack  model  (figure  3.12)  gives  reasonable  variance  reduction  of 
about  41%.  Figure  3.12a  shows  that  small  misfits  are  achieved  for  symmetry 
axis  locations  in  two  regions  in  the  Southwest  and  Northeast.  Only  the 
Southwestern  region  gives  rise  to  statistically  significant  solutions.  The  lowest 
contom  shows  the  confidence  region  where  the  critical  misfit  level  is  computed 
from  the  f-distribution  as 

ScH*  =  So,t[l  +  J^fk,N-k{l-a)]  (25) 
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Free-Surface  Effect  on  Polarization  IncAngle 


Figure  3.4:  Free-sudace  effect  on  pdamation  incidence  angles  for  a  set  of 
Poisson  ratios. 
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2.  MONITORING  NUCLEAR  TEST  SITES 
WITH  GERESS 


M.  L.  lost,  J.  Schweitzer,  and  H.*P.  Harjes 


2.0  ABSTRACT 


The  GERESS  airay  in  south-eastern  Germany  has  been  proposed  to  parti¬ 
cipate  as  an  a-stadon  during  the  upcoming  GSETT-3  test  With  the  aim  of 
monitoring  a  comprehensive  test  ban  treaty,  special  interest  has  emerged  to 
asses  the  capability  of  this  new  array  to  monitor  the  known  nuclear  test  sites. 
During  the  GERESS  site  survey  (1988-1989)  and  4  1/2  years  of  GERESS 
operation  (Jan.  1990  -  June  1S>94),  44  presumed  nuclear  events  have  been 
recorded:  18  from  Nevada,  14  from  the  Tuamotu  Archipelago,  6  from  eastern 
Kazakhstan,  5  from  Lop  Nor,  and  1  from  Novaya  Zemlya.  The  detection  thres¬ 
holds  at  GERESS  are  01^  =  4.3  for  Nevada,  mb  =  4.0  for  Tuamotu  Archipelago, 
mb  =  4.6  for  Lop  Nor,  nib  *=  4.3  for  Novaya  Zemlya,  and  tUb  =  4.S  for  eastern 
Kazakhstan.  For  Nevada  and  the  Tuamotu  Archipelago,  these  detection  thres¬ 
holds  at  GERESS  cotiespond  to  a  yield  estimate  in  the  order  of  1  kt  at  each 
site.  These  path  specific  detection  thresholds  deviate  substantially  from  those 
derived  from  global  anqilitude-distance  relations  (i.e.,  Gutenberg-Richter). 
Based  (»  a  cross-conelation  method  and  a  cluster  analysis,  GERESS  records 
can  be  used  to  identify  specific  test  sites. 


2.1  INTRODUCTION 

The  Group  of  Scientific  Experts  (GSE)  of  the  Conference  of  Disarmament 
(CD)  in  Geneva  has  started  d^eloping  a  final  concept  for  an  international 
seismic  data  exchange  system  (CD/1211,  1993)  providing  technical  support  for 
monitaring  a  coinprehensive  nuclear  test  ban  treaty  (CTBT).  The  near  real 
time  system  includes  a  three-tiered  global  network  of  seismic  staticms  where  the 
first  tier,  the  alpha-network,  consists  about  half  of  arrays  and  half  of  three- 
component  stations.  These  a^ha-stations  are  designed  to  provide  not  only  the 
leqitired  detection  threshold  throu^iout  die  world  but  also  to  achieve  a  prelim¬ 
inary  event  location  s^ch  is  then  inqaoved  by  the  supplementary  data  of  the 
beta  and  gamma  networks.  In  diis  respect,  the  a-netwoik  of  about  60  globally 
distriboted  stations  is  regarded  as  a  teleseismic  monitoring  system. 
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Table  2-1:  Distances  of  Nuclear  Test  Sites  to  GERESS 


Country 

Test  Site 

Lat 

Long 

A 

[deg] 

BAZ 

[deg] 

China 

Lop  Nor 

41.6  N 

88.7  E 

51.3 

68 

Fiance 

Fangataufa 

222  S 

138.8  W 

145.6 

311 

France 

Mururoa 

21.9  S 

139.0  W 

145.4 

311 

Great  Britain 

Nevada 

37.2  N 

116.2  W 

83.5 

322 

USA 

Nevada 

37.2  N 

116.2  W 

83.5 

322 

USSR 

Kazakhstan 

49.9  N 

78.8  E 

41.2 

63 

USSR 

Novaya  Zemlya 

73.4  N 

54.7  E 

30.4 

22 

39 


Figure  2-5:  Same  as  Fig.  2-4,  except  filter:  0.7  -  2.0  Hz. 
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3.  P-WAVE  ARRAY  POLARIZATION  ANALYSIS  AND 
EFFECTIVE  ANISOTROPY  OP  THE  BRITTLE  CRUST* 

G.  H.R.  Bokelmann 


3.0  ABSTRACT 


P-wave  polarization  constrains  local  anisotropy  in  the  vicinity  of  the  re¬ 
ceivers.  Using  three-component  and  array  data  from  the  regional  GERESS 
zirray  in  Southeastern  Germany,  we  measme  polzuization  a  and  propaga¬ 
tion  vectors  s  for  P-phases  of  120  events.  Angular  deviations  a-s  between 
these  normalized  vectors  often  approach  10“,  rendering  them  easily  measur¬ 
able.  The  effect  of  anisotropy  can  be  distinguished  &om  remote  effects,  since 
all  remote  effects  such  as  source  mislocation,  distant  lateral  heterogeneity 
or  distant  anisotropy  affect  polarization  and  propagation  vectors  simultzine- 
ously.  Averaging  removes  sensitivity  to  near-receiver  heterogeneity,  and  local 
Einisotropy  is  left  as  sole  cause  of  the  effects  in  a-s.  This  method  hence  gives 
local  effective  anisotropy  in  the  near-receiver  cnxst  averaging  over  a  depth 
interval  of  a  wavelength  (ss  6km).  We  resolve  strike  and  dip  of  the  symetry 
plane  and  also  two  dimensionless  numbers  t}  and  r  which  give  constraints  on 
four  of  the  elastic  parameters.  The  optimum  model  (variance  reduction  44%) 
has  symetry  plane  orientation  of  strike  113*  and  dip  49*  to  the  North,  which 
corresponds  closely  to  the  consistently  observed  gneiss  foliation  direction  in 
the  area  (120*,  50*-60*).  Comparing  t)  and  r  with  predictions  &om  different 
physical  models  we  find  that  the  data  are  fit  by  a  gneiss  model  assuming  that 
the  anisotropy  is  dominated  by  the  mica,  if  3%  to  8%  of  the  mica  are  well- 
aligned.  This  suggests  that  anisotropy  in  the  region  studied  is  dominated  by 
the  effect  of  local  foliation  rather  than  the  regional  stress  field. 

Key  words:  array  seismology,  wave  propagation,  polarization,  anisotropy, 
inversion 

3.1  INTRODUCTION 


In  the  pioneering  study  of  anisotropy  in  the  Earth  Hess  (1964)  showed 
that  azimuthal  travel  time  delays  of  oceanic  P„  can  be  explained  by  aligned 
olivine  crystals  in  the  oceanic  upper  mantle.  For  such  a  phase,  which  travels 
subhorizontally  on  most  of  its  path.  Backus  (1965)  showed  that  azimuthal 
travel  time  variation  is  characterized  by  a  rather  simple  trigonometric  ex¬ 
pression,  given  that  anisotropy  is  uwak.  Fitting  of  such  trigonometric  terms 
has  subsequently  been  used  for  teleseismic  P-waves,  f.e.  in  Dziewonski  and 
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(Jenkins  and  Watts,  1969),  where  Sopt  is  the  minimum  misfit  and  k=4  the 
number  of  parameters.  Higher  contours  give  2  and  4  times  the  separation  of 
the  critical  level  (1  —  a  =  0.95%  confidence  level)  firom  the  minim^  level. 
The  optimum  model  has  a  symmetry  axis  orientation  of  azimuth  6  =  208° 
and  dip  $  =  36°.  The  symmetry  plane  is  oriented  normal  to  the  symmetry 
ends,  and  has  strike  118°  and  dip  54°  to  the  North.  The  optimxim  model 
corresponding  to  the  global  minimum  in  the  Southwest  is  shown  in  figure 
3.12b  by  its  velocity  and  polarization  pattern.  To  be  directly  comparable  with 
the  data,  only  the  partial  lower  hemisphere  from  0°to  55°is  shown.  Beside  the 
symmetry  plane,  also  the  symmetry  axis  shows  up  with  fast  velocity  in  the 
lower  part  of  the  figure.  The  polarization  deviations  on  the  right-hand  side 
of  figure  3.12b  are  given  for  0°,  10°,  ...,  50°incidence  and  may  be  comp2u:ed 
with  the  data  in  figtire  3.7.  The  major  features  of  the  data  set  are  explained 
by  the  model. 

The  results  for  assumed  dry  and  sattirated-crack  model  are  quite  different. 
The  saturated-crack  model  gave  reasonable  fit  to  the  data  and  a  strike  angle 
of  the  high-velocity  plane  of  118°,  while  the  dry-crack  was  in  clear  confiict 
with  the  data. 

Now  we  perform  the  more  general  nonlinear  search  for  the  parameter  set 
{6,  $,  Tj,  t),  which  unlike  the  previous  example  does  not  assume  a  pEuiicular 
physical  fracture  model.  With  4  parameters  and  the  given  botmds  on  rj  and  r 
this  search  is  not  particularly  compute-intensive,  but  guarantees  the  globally 
optimal  solution,  given  that  the  bounds  are  appropriate  and  the  sampling 
is  sufficiently  fine.  The  misfit  pattern  in  figure  3.13a  is  not  unlike  the  one 
for  the  satxirated  crack  model.  However,  the  optimum  in  the  Southwest  is 
more  pronounced  now  and  the  variance  reduction  is  three  percent  larger, 
now  about  44%.  The  resulting  parameters  are  0  =  203.3°,  $  =  41.05°, 
T]  =  0.887  and  r  =  1.09.  The  maximum  P-velodty  difierence  of  the  model 
is  about  4.7%.  The  change  in  orientation  from  the  fit  of  model  HCSl  by  a 
few  degrees  is  insignificant  considering  the  large  confidence  region  in  figme 
3.12a.  In  the  general  model  the  confidence  region  is  substantially  smaller. 
For  the  polarizations,  the  main  difference  is  the  behaviour  in  the  South- 
Western  part,  where  the  general  model  predicts  a  somewhat  smaller  effect. 
This  corresponds  to  the  smaller  velocity  contrast  in  that  region.  The  high- 
velocity  symetry  plane  strikes  with  about  113°  and  dips  with  about  49°  to 
the  Northeast. 

In  the  search  we  varied  ri  and  r  independently.  But  can  they  be  resolved 
independently  by  the  data  constraints  or  are  they  subject  to  a  major  tradeoff? 
Figure  3.14  shows  the  misfit  aroimd  the  optimum  model  keeping  the  symetry 
axis  orientation  fixed  and  varying  rj  and  r.  The  lowest  contour  again  shows 
the  confidence  region  of  admissible  tj  and  r,  which  extends  for  r)  from  0.76 
to  1.07  and  for  r  from  1.  to  1.2  .  The  isotropic  case  at  r;  =  1  and  r  =  1  is 
well  outside  the  highest  contour  for  four  times  the  critical  misfit  separation. 

The  strong  ascent  for  A'  near  0  strongly  discriminates  against  the  dry 
crack  model. 
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Recently,  the  GERman  Experimental  Seismic  System  (GERESS)  has  been 
pnqx>sed  to  participate  in  the  upcoming  GSE  Tech^cal  Test  3  (GSETT-3)  as 
an  a-station.  Hence,  the  ciq)abilities  of  the  GERESS  array  to  monitor  the 
presently  known  nuclear  test  sites  has  become  of  interest 

Based  on  4  1/2  years  of  GERESS  data  (Jan.  1990  -  June  1994),  this  piq)er 
addresses  the  detection  «q>abilities  fm  presumed  nuclear  events.  We  also  ad^ 
10  events  from  the  GERESS  site  surveying  phase  (1988-1989)  to  include  data 
from  the  joint  verification  experiments  (JVE)  and  from  the  Kazakhstan  test  site. 
Due  to  the  amount  of  data,  we  focus  on  events  at  the  Nevada  Test  site  (NTS) 
and  events  in  the  Tuamotu  Archipelago. 


2.2  THE  GERESS  ARRAY 


Since  the  establishment  of  the  Giafenberg  array  (GRI^  in  1975,  consider¬ 
able  efforts  have  been  made  in  Germany  to  use  array  data  far  mmutoring 
underground  nuclear  explosions  (Haijes  and  Hanka.  1986).  After  it  became 
clear  that  high  frequency  signals  play  an  inqxrrtant  role  for  detection  and 
identification  of  small  underground  nuclear  tests  at  regional  distances  (e.g., 
Mykkeltveit  et  al.,  1990;  Ringdal,  1990),  a  2S-element  shon  period,  small  aper¬ 
ture  array  in  south-eastern  Germany  was  planned.  This  array  c(xiq>lcments  the 
broad-band,  large  aperture  Grafenberg  array  which  had  not  been  designed  for 
detection  purposes.  The  design  goals  of  GERESS  included  the  lowering  of  the 
global  detectirm  threshold  in  addition  to  nx>nitc»ing  local  and  regional  seismi¬ 
city  in  Central  Europe  and  the  Mediterranean  earthquake  zone. 

The  siting  survey  for  diis  new  array  began  in  1987  (Haijes,  1989).  During 
this  survey,  single  stations  were  set  up  near  die  localities  Gross-Lichtenberg 
(GRLB  48.838  N,  13.720  E),  Haidmuehle  (HAID  48.800  N,  13.786  E),  and 
Sulzberg  (SULZ)  at  the  present  GERESS  C2  location.  These  single  stations 
operated  event  triggered  from  April-September  1988.  From  December  1988  to 
/qitil  1989,  a  9-element  test-array  was  in  cqieration  at  the  present  GERESS  site 
(Hatjes,  1990a).  The  25-eiement  GERESS  array  was  installed  in  cooperati<» 
with  Southern  Methodist  University  (Dallas,  Texas)  and  Ruhr-University 
Bochum  and  records  are  available  from  Apr.  1990  (Haijes,  1990b). 

GERESS  is  located  in  the  Bohemian  Forest  in  south-eastern  Germany  near 
the  border  to  Austria  and  the  Czech  Republic  (Hg.  2-1).  The  array  consists  of 
25  stations  with  vertical-coriqionent  short-period  Teledyne  Geotech  GS-13 
instruments  sairqiled  at  40  Hz.  In  addidrm,  four  of  the  sites  include  horizontal 
component  instruments.  At  the  key  station  of  die  array,  GEC2  (48.84511  N, 
13.70156  E,  1132  m),  a  STS-2  bioad-band  instrument  is  saiiqiM  at  80  Hz. 
The  geometry  of  die  array  is  based  <»  concentric  rings  proving  an  overall 
qieituie  of  about  4  km  (Fig.  2-2).  Data  from  the  array  are  continuously 
transmitted  to  the  NORSAR  data  center  in  Kjeller  and  to  tte  Institute  of  Geo¬ 
physics,  Ruhr-University  Bochum. 
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Tabic  2-2:  Parameters  of  Events  from  Nevada  Test  Site  Recorded  at  GERESS  1988-1993 
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4:  Amplitude  measured  on  trace  filtered  1.0- 3.0  Hz 


Anderson  (1983).  Interpretation  of  such  data,  however,  is  complicated  by  the 
fact  that  lateral  heterogeneity  can  also  cause  azimuthal  variation.  This  trade¬ 
off  is  often  difficult  to  resolve  from  travel  time  data  alone.  The  second  type 
of  "classical”  sinisotropy  studies  is  shear- wave  splitting  (Ando,  Ishikawa  and 
Yamazaki,  1983;  Silver  and  Chan,  1991),  a  technique  which  uses  the  strong 
effect  of  zmisotropy  on  waveforms  of  shear-waves  due  to  different  speed  of 
the  two  orthogonally  polarized  shear-waves.  The  constraint  on  anisotropy  is 
not  unlike  that  of  travel  times.  Both  give  integral  constraints  only,  namely 
on  bulk  anisotropy  along  the  wave  path.  Hence,  for  both  methods  a  tradeoff 
exists  between  size  of  anisotropy  and  extent  of  anisotropic  path.  More  im¬ 
portantly,  an  anisotropic  region  is  difficult  to  localize  using  these  data.  For 
example,  msintle  anisotropy  csm  be  inferred  ffom  shear-wave  splitting  data 
only  if  observed  effects  are  too  large  to  be  explained  by  crustal  anisotropy 
alone.  This  is  somewhat  unsatisfying,  and  methods  for  independent  deter¬ 
mination  of  anisotropy  in  the  crust  tmder  the  receivers  should  be  sought. 
This  paper,  in  fact,  demonstrates  that  beside  active  experiments  an  entirely 
independent  constraint  on  near-receiver  anisotropy  is  possible;  in  this  case  it 
follows  from  analysis  of  three-component  array  data. 

Anisotropy  studies  using  travel  times  and  shear-wave  splitting  can  be  per¬ 
formed  with  single-station  data,  in  the  later  case  requiring  a  three-component 
instrument.  For  the  present  study  of  P-wave  phases,  we  use  array  and  three- 
component  data  from  the  GERESS  array.  These  array  data  allow  determi¬ 
nation  of  the  propagation  vector  s  =  k/u>  of  plane  waves  crossing  the  array. 
While  this  particular  feature  of  beamforming  is  well-appreciated  for  signal 
enhancement,  event  detection  and  localization  (Beauchamp,  1975),  this  par¬ 
ticular  strength  has  apparently  not  been  used  in  anisotropy  studies  so  far. 

The  propagation  vector,  given  by  s  =  Vt,  may  be  specified  by  azimuth  S, 
and  horizontal  slowness  p  =  -y/s**  -I-  or  alternatively  incidence  angle  = 
sin~^(jn})  with  the  local  velocity  v.  On  the  other  hemd,  three-component 
stations  adlow  determination  of  the  polarization  vector  a,  which  gives  the 
dominant  direction  of  linear  particle  motion,  a  can  be  specified  by  azimuth 
6a  and  incidence  angle  <f>a.  In  the  following  we  restrict  our  intention  to 
normalized  s  and  a,  namely  to  azimuth  and  incidence  angles.  This  effectively 
ignores  that  the  zuray  also  gives  estimates  of  the  true  or  apparent  velocity. 
Neither  of  the  two  azimuth  estimates  9,  and  9a  necessarily  correspond  to  the 
azimuth  of  the  receiver-source  great  drde,  since  heterogeneity,  anisotropy 
and  other  effects  may  cause  deviations.  The  same  is  true  for  the  incidence 
angles.  For  known  source  location  we  get  predictions  (9f,  <f>p)  for  a  reference 
Earth  model.  Comparing  (9,,^,)  and  (9a,  ^a)  with  the  predictions  gives  rise 
to  so-called  mislocation  diagrams.  For  the  purpose  of  event  localization  the 
two  data  types  and  correspondingly  their  mislocation  patterns  are  typically 
treated  the  same,  where  the  azimuth  is  used  for  localization  piuposes,  while 
the  incidence  angle  serves  for  identification  of  the  phase  for  regional  events 
and  epicentral  distance  in  the  case  of  teleseismic  events. 

Interestingly,  source  location  procedures  have  so  far  disregarded  that  these 
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3.6  WHAT  CAUSES  THE  EFFECTIVE  ANISOTROPY? 

In  principle,  intrinsic  anisotropy  as  well  as  macroscopic  effects  from  ori¬ 
ented  cracks,  fractures  or  layering  rjtn  catise  anisotropy  of  this  size.  To 
understand  the  nature  of  this  anisotropy,  it  is  of  large  value  to  consider  ge¬ 
ological  evidence  for  the  region.  The  array  is  located  on  outcrop  of  granite 
and  gneiss.  While  the  granite  apparently  has  no  preferred  alignment,  the 
gneiss  shows  strong  foliation  with  consistent  orientations.  The  strike  is  ap¬ 
proximately  120‘’and  the  dip  50-60‘’(0tt,  pers.comm.).  This  closely  coincides 
with  the  orientation  of  the  high-velocity  planes  for  the  best-fitting  model  from 
the  imconstrained  inversion  (figure  3.13)  suggesting  that  rock  foliation  is  the 
major  cause  of  anisotropy  in  this  case.  The  size  of  the  anisotropy  suggests 
that  the  foliation  (hrection  is  spatially  consistent  over  some  distance,  perhaps 
several  kilometers.  In  fact,  surface  geology  does  suggest  that  this  is  indeed 
the  case  (Ott,  pers.comm.).  Rock  foliation  as  a  cause  of  effective  anisotropy 
has  been  discussed  before  (f.e.  Luschen  et  al.,  1991). 

But  is  rock  foliation  the  only  factor  involved  and  how  does  the  overall 
anisotropy  relate  to  intrinsic  properties  of  minerals?  Gneisses  typically  con¬ 
tain  large  amounts  of  mica,  a  mineral  giving  rise  to  very  large  anisotropy. 
Since  these  mineral  are  usually  quite  well-oriented  (Wang  et  al.,  1975),  we 
may  expect  major  effect  on  the  observed  anisotropy.  In  fact,  mica  orienta¬ 
tions  are  typically  dominant  in  forming  macroscopic  gneiss  foliation  planes 
(Suppe,  1985).  Pure  mica,  however,  does  not  explain  the  observed  effect. 
Aleksandrov  and  Ryzhova  (1961)  give  elastic  constants  for  muscovite,  for 
which  we  compute  r}  =  0.625  and  r  =  3.88  .  The  latter  value  is  clearly 
incompatible  with  our  results.  More  realistic  modds  would  include  struc¬ 
tural  complexity,  the  occurence  of  several  types  of  intrinsic  anisotropy,  and 
fractures  or  cracks  in  the  subsurface.  In  fact,  we  have  physical  models  at 
our  disposition  predicting  some  of  these  effects.  For  a  medium  consisting 
of  anisotropic  layers,  Schoenberg  and  Douma  (1988)  give  expressions  for  av¬ 
eraging  elastic  constants  to  obtain  effective  elastic  constants  in  the  long- 
wavelength  assumption.  For  such  a  medium,  also  effective  parameters  fj  and 
f  are  composed  of  the  averaged  elastic  constants  as 
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summarizes  the  GERESS  recordings  of  underground  nuclear  explosions.  Path-specific  delectioo  diiesholds  at  GERESS  are 
derived  for  NTS,  Tuainolu  Archipelago,  Kazakhstan,  Novaya  Zeinlya,  and  Lop  Nor.  These  empirical  detection  thresholds 
substantially  deviate  from  those  raimiatwi  from  global  amplitude-distance  relations  (Le.  Gutenbog-Richter).  This  case  history 
again  casts  some  doubts  on  theoretical  network  capability  estifnates  unless  these  are  based  on  careful  station  calibration. 

The  second  paper  (Btdtelmann)  uses  three-component  and  conventional  vertical-component  array  data  from  GERESS  to 
separately  measure  polarization  and  slowness  vectors.  This  data  set  allows  to  distinguish  the  effect  of  anisotropy  from  lateral 
heterogeneity,  from  interestiiig  geophysical  conclusions  about  the  origin  of  this  anisotn^y,  the  study  has  also  implica- 
tions  on  the  combined  use  of  directknal  data  from  arrays  and  three-component  stations  in  a  global  monitoring  netwoik.  As 
local  anisotropy  in  the  vicinity  of  the  receivers  influences  P-wave  polarization,  there  may  result  a  systematic  bias  between 
slowness  parameters,  measured  with  arrays  and  direcrional  data,  derived  from  tlnee-component  statkws. 
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The  GERESS  data  center  at  Bochum  is  operating  an  experimental  on-line 
processing  system  (lost  1992b),  which  uses  software  devel(q>ed  at  NORSAR 
(Mykkeltveit  and  Bungum,  19^;  Fyen,  1987)  and  modified  in  Bochum.  The 
autmnatic  system  detects  teleseismic  phases  in  additiem  to  locating  regional 
events  (Haijes  et  al.,  1993;  lost,  1993).  Recently,  Haijes  et  al.  (1994) 
evaluated  the  contribution  of  GERESS  to  global  seismicity  using  the  NEIC  bul¬ 
letin  as  reference  and  concluded  that  GERESS  has  excellent  teleseismic  detec¬ 
tion  a^rabilities. 

GERESS  is  located  in  the  teleseismic  distance  range  with  respect  to  all 
nuclear  test  sites  that  were  active  1988  -  June  1994  (Fig.  2-3).  Table  2-1  indi¬ 
cates  the  mean  locations  of  these  test  sites,  the  qricentral  distances  to  GERESS, 
and  the  back-azimuths.  This  distance  range  may  actually  be  no  drawback  for 
the  task  of  mcmiKning  a  low-yield  test  ban  treaty.  For  example,  GERESS 
autmnatically  detected  a  P-wave  signal  from  a  chemical  calibration  experiment 
of  1  kt  charge  at  NTS  (smnetimes  called  non-proliferation  test  (NFT))  whereas 
Lg  waves  were  barely  visible  at  the  low  noise  statimi  LTX  at  regional  distance 
(13*’,  Lawson,  pers.  comm.). 


NEVADA  TEST  SITE 


Frexn  the  total  of  24  underground  nuclear  explosions  conducted  during  the 
investigated  time  period,  18  wavefcHins  are  available  at  GERESS  in  addition  to 
the  aforementioned  chemical  NPT.  4  of  the  24  events  (1990^/20  9:15  local, 
1990/9/27  11:02  local,  1992/4/30  9:30  local,  and  1992/6/19  VICTORIA) 
were  too  small  to  {uoduce  any  observable  ngnal  at  stations  in  Ontral  Europe. 
During  2  tests  (METROPOLIS  aixl  DISTANT  ZENITH),  GERESS  was  out 
Table  2-2  shows  parameters  of  events  at  NTS  from  1988-1993.  The  two  1988 
events  COMSIXXIIK  and  KEARSARGE  (JVE)  were  teemded  dining  the  site 
survey  of  GERESS  at  station  GRLB. 

Figure  2-4  shows  the  waveforms  of  all  observed  events  from  Nevada, 
sorted  according  to  magnitndf.  determined  at  GERESS  (Table  2-2).  All  traces 
start  10  seconds  before  the  theoretical  onset  time  calculated  from  the  published 
Department  of  Energy  (DOE)  parameters  as  given  in  NEIC  using  the  IASP91 
(Kennett  and  Engdahl,  1991)  travel  time  tables.  For  the  2  events  in  1988,  sin¬ 
gle  traces  are  shown.  For  the  evoits  in  1990-1993,  the  qrtimum  GERESS  beam 
is  displayed.  In  Hgures  2-4  -  2-6,  we  filtered  all  available  vertical  traces  with  a 
third  order  butterwortfa  band-pass  filter  from  0.5  Hz  -  2.5  Hz,  0.7  Hz  -  2.0  Hz, 
and  1.0  Hz  -  3.0  Hz,  respectively,  before  forming  die  displayed  optimum  beam. 
Only  forward  filtering  is  used  throughout  this  paper.  The  parameters  (back- 
azimuth  and  slowness)  to  calculate  the  tqitimum  GERESS  beam  were  obtained 
from  fk  arudysis  using  all  available  vertied  traces  at  GERESS. 

In  Figures  2-4  -  2-6,  we  clearly  see  a  scatter  of  +/-  1  sec  for  die  observed 
onset  times  at  GERESS.  A  first  eiqilanation  would  be  duu  the  abstdute  timing 
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Table  2-3:  Magnitude  Yield  Relations  for  NTS 


Yr 

Mo 

Day 

Name  /  Site 

“b 

GRF 

Yield 

GRF 

“b 

GERESS 

Yield 

GERESS 

Yield 

Nuttli 

No. 

1988 

Jun 

2 

Comstock 

5.4 

60 

5.4 

49 

68 

13 

1988 

Aug 

17 

Kearsarge,  JVE 

5.5 

75 

5.4 

150 

68 

12 

1990 

Jun 

13 

Bullion 

5.8 

150 

5.7 

114 

150 

19 

1990 

Jun 

21 

Austin 

<3.9 

<0.8 

2 

1990 

Jul 

25 

Mineral  Quarry 

5.0 

20 

4.7 

7 

8 

9 

1990 

Oct 

12 

Tenabo 

5.6 

90 

5.5 

65 

92 

16 

1990 

Nov 

14 

Houston 

5.5 

70 

5.4 

49 

68 

14 

1991 

Mar 

08 

Coso 

4.2 

2 

7 

1991 

Apr 

04 

Bexar 

5.7 

no 

5.5 

65 

92 

15 

1991 

Apr 

16 

Montello 

5.4 

50 

5.3 

37 

50 

10 

1991 

Aug 

15 

Floydata 

<3.6 

<0.3 

1 

1991 

Sep 

14 

Hoya 

5.6 

90 

5.6 

86 

125 

18 

1991 

Oct 

18 

Lubbock 

5.2 

30 

5.4 

49 

68 

11 

1991 

Nov 

26 

Bristol 

4.9 

10 

4.7 

7 

8 

8 

1992 

Mar 

26 

Junction 

5.6 

100 

5.6 

86 

125 

17 

1992 

Jun 

23 

Galena 

<4.0 

<1 

4 

1992 

Sep 

18 

Hunters  Trophy 

4.2 

2 

6 

1992 

Sep 

23 

Divider 

<3.9 

<0.7 

3 

1993 

Sep 

22 

NPT 

4.0 

1 

5 

The  GERESS  yields  are  from  log  y  «  1.21  mb  -  4.84.  The  Nuttli  yields  were  calculated  using 
his  1986  paper  and  GERESS  magnitudes  (see  text). 

The  No.  column  indicates  a  sotting  of  the  events  according  to  mb(GERESS). 
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Figure  2-7.  Correiudra  matrix  of  the  events  from  Nevada  test  site  (for  event 
numbers,  sec  Table  2-3). 
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EVENT  NUMBER  (FLOYDATA  =  1.  BULLION  =  19) 
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two  observed  quantities,  namely  propagation  and  polarization  vectors,  are 
principally  different.  While  they  coincide  in  an  isotropic  Earth,  elastic 
anisotropy  in  the  Earth  would  affect  propagation  stnd  polarization  vectors 
differently.  For  event  localization,  this  means  that  there  is  a  systematic  bias 
between  propagation  and  polarization  vectors,  which  should  be  removed  by 
taking  local  anisotropy  into  account.  This  systematic  bias,  however,  is  the 
prime  quantity  of  interest  in  this  study,  since  it  offers  us  a  new  constraint 
on  anisotropy  in  the  vicinity  of  the  recovers  (in  fact  exclusively  there).  In 
this  study  we  use  polarization  and  propagation  vectors  exclusively.  Travel 
time  and  amplitude  information  is  not  used,  since  that  might  obscure  the 
near-receiver  crustal  anisotropy  with  information  gathered  somewhere  along 
the  ray  paths. 


3.2  OBSERVATIONS 


3.2.1  GERESS  Array 

The  GERman  Experimental  Seismic  System  GERESS  (Harjes,  1990), 
shown  in  figure  3.1,  is  a  regional  seismological  array  with  aperture  of  about 
4  km  consisting  of  25  stations  with  1  Hz  vertical  short-period  instruments 
(Geotech  GS-13)  sampled  at  40  Hz,  with  4  of  them  also  equipped  with  hor¬ 
izontal  components  (A2,D1,D4,D7;  triangles).  Additional  three-component 
instruments  at  the  C2  location  are  sampled  at  120  Hz  (Geotech  GS-13)  and 
10  Hz  (Geotech  BB-13).  For  further  descriptions  of  installation  criteria,  per¬ 
formance  and  operation  see  Harjes  (1990)  and  Harjes  et  al.  (1993).  The 
array  is  located  in  Eastern  Bavaria,  Germany,  in  a  hilly  region  on  crystalline 
outcrop  of  the  Bohemian  massive.  The  topographic  variation  across  the  ar¬ 
ray  is  about  200  meters.  A  major  advantage  of  data  &om  the  GERESS  array 
is  that  there  are  no  sediments  under  the  stations,  greatly  improving  the  qual¬ 
ity  particularly  of  the  polarization  data.  Otherwise  P-S  conversions  at  the 
crust-sediment  boimdary  might  disturb  polarizations  substantially. 

3.2.2  Events 

This  study  uses  a  set  of  95  events  recorded  at  GERESS  array  stations  in 
the  time  interval  between  May  1990  and  January  1993  (table  3.1).  These 
events  comprise  the  whole  range  &oxn  re^onal  to  telesrismic  distance.  The 
emphasis  in  this  study  is  on  first-arriving  phases.  Only  the  later-arriving 
phase  PcP  is  included  in  three  cases.  Ihble  3.1  is  part  of  a  larger  set  of 
events.  For  the  events  not  included  in  the  table  either  the  polarization  or 
propagation  vector  had  to  be  discarded,  which  was  the  case  particularly 
often  early  in  the  experiment.  In  general,  events  were  discarded  if  data 
&om  less  than  3  three-component  stations  or  9  array  stations  were  available. 
Out  of  112  events  with  acceptable  polarization  data,  81  had  data  for  all  6 
three-component  instruments  present.  Out  of  108  events  with  acceptable 
propagation  data,  88  had  data  from  at  least  20  array  stations. 

Next  we  discuss  procedures  of  extracting  backazimuth  and  incidence  angle 


C44  =  1/  <  1/C44  > 


(31) 


where  <>  denotes  the  thickness- wdghted  average  of  the  respective  quantity. 
If  we  take  the  simple  approach  of  assuming  that  due  to  their  large  effect  only 
mica  contribute  to  the  bulk  anisotropy,  we  may  simply  compute  rj  and  r  for 
different  percentages  of  mica  content.  For  a  mica  fraction  of  30%  mica  after 
Wang  et  al.  (1975)  and  70%  other  material  assumed  to  be  isotropic  we  obtain 
T)  =  0.816  and  r  =  1.74  .  While  this  is  far  outside  the  confidence  region  of 
the  observation,  values  for  3%  to  8%  mica  content  fall  within  the  confidence 
region  (figure  3.14).  These  are  low  values  compared  with  mica  contents  in 
typical  gneisses  (Wang  et  al.,  1975).  We  have  to  keep  in  mind  though  that 
certainly  not  all  of  the  mica  minerals  are  fuUy  aligned.  Other  constituents 
may  also  play  a  minor  role  in  defining  the  bulk  anisotropy. 

Satisfying  the  polarization  data,  this  gndss  model  is  also  a  plausible  ex¬ 
planation  in  conjtmction  with  the  geological  data.  Another  possible  model, 
however,  is  that  of  cracks  or  fractures  aligned  in  the  direction  of  planes  of 
weakness  which  we  coincide  with  the  gneiss  foliation.  Fitting  a  model  of 
aligned  cracks  after  Hudson  (1981)  had  given  essentially  the  szune  orienta¬ 
tion.  Resulting  parameters  were  strike  114®  and  dip  50°  for  the  symetry 
planes.  We  also  obtained  a  crack  density  of  e  =  0.09  and  two  Lame  parame¬ 
ters  of  the  crack  material  A'  =  0.049  GPa  and  /i'  =  0.0  GPa.  The  associated 
error  bars  discriminate  clearly  against  the  dry-crack  case,  where  A'  =  O.GPa 
.  Such  dry  cracks  give  rise  to  large  values  of  r,  which  we  have  previously 
argued,  are  in  conflict  with  the  data.  Based  on  the  achieved  misfit,  the  model 
of  cracks  opening  along  the  foliation  planes  can  not  be  niled  out. 

It  appears  difficult  to  distinguish  between  models  of  intrinsic  anisotropy 
and  cracks  or  fractures.  In  this  case,  one  should  consider  all  possible  models. 
In  fact,  the  occurence  of  cracks  along  planes  of  weakness  given  by  foliation 
planes  is  not  iinlikely.  For  small  pressures  in  the  laboratory,  it  is  very  com¬ 
mon  to  observe  opening  of  cracks  along  the  foliation  direction.  This  can  be 
demonstrated  f.e.  with  the  data  of  Wang  et  al.  (1975).  From  their  data  we 
obtain  values  of  r  »  1.25  at  Ikbar  pressure  and  much  larger  values  up  to 
10  at  lower  pressure.  Opening  of  cracks  can  be  caused  by  the  local  stress 
field.  Values  for  the  strike  of  the  maximum  horizontal  stress  direction  an  in 
the  area  are  around  149°  in  the  KTB  borehole  (Brudy,  Fuchs  and  Zobeurk, 
1993),  but  there  are  also  values  in  the  range  of  115°  (f.e.  Falkenberg  granite, 
Rummel  et  al.,  1983).  These  values  are  not  far  from  the  foliation  strike. 
Hence  crack  opening  along  the  foliation  direction  is  not  unlikely,  particularly 
if  the  fluid  pressure  were  high.  However,  we  note  that  anisotropy  in  our  data 
set  appears  to  be  due  to  foliation  rather  than  vertical  cracks  perpendicular 
to  the  smallest  horizontal  stress  direction.  On  the  other  hand,  the  regional 
stress  field  may  have  an  effect  due  to  opening  of  cracks  preferentially  parallel 
to  the  planes  of  weakness,  which  are  given  by  the  gneiss  foliation  planes. 

Perhaps,  more  insight  can  be  gained  if  also  shear  waves  are  used.  Figure 
3.15  gives  the  prediction  of  the  shear  wave  velocities  for  the  optimum  model. 
However,  unambiguous  determination  of  upper  cnistal  anisotropy  requires  in- 
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of  die  GERESS  array  was  not  stable  and  imxiuced  these  deviations.  But  this 
c(»clusion  could  be  rejected  because  all  GERESS  onset  times  show  a  constant 
travel  time  difference  of  7.S  -  7.8  sec  with  respect  to  published  onset  times  of 
these  events  at  the  Grafenberg  airay.  It  is  obvious  that  the  same  absolute  time 
scatter  could  not  occur  at  two  different  arrays.  In  addition,  we  verified  the  rela¬ 
tive  (mset  times  of  MONTELLO  and  HOUSTON  at  stations  Berggiesshiibel 
(BRG),  Colltnbeig  (CLL),  and  Kaspersk’e  Wary  (KHC).  The  observed  scatter 
of  (mset  times  can  be  caused  in  part  by  the  velocity  structure  at  NTS 
(McLaughlin  et  al.,  1987;  Ferguson  et  aL,  1989).  Because  a  time  dependent 
travel  time  between  NTS  and  the  two  arrays  in  Germany  is  not  realistic  our 
preferred  hypothesis  is  that  the  published  source  times  are  not  correct 

For  the  events  in  1990-1993,  we  measured  arrival-times,  amplitudes,  signal 
to  noise  ratios  (SNR),  and  periods  on  the  beam  traces  (Table  2-2).  It  is  not 
surprising,  diat  very  small  slowness  values  (the  theoretical  value  finn  IASP91 
is  5.13  sec/deg)  were  determined,  since  the  GERESS  aperture  of  about  4  km 
limits  the  resolution  of  slowness  measurements  in  the  teleseismic  distance 
range. 

The  nuclear  tests  AUSTIN,  FLOYDATA,  GALENA,  and  DIVIDER  were 
not  detected  by  the  automatic  STA/LTA  detector  at  GERESS  but  show  weak 
signals  at  their  expected  onset  times  in  Figures  2-4  -  2-6.  The  measured  param¬ 
eters  (Table  2-2)  indicate  an  upper  limit  for  tiiese  events.  In  the  filter-bt^  0.7 
-  2.0  Hz  (Fig.  2-5),  weak  cmsets  of  AUSTIN,  GALENA,  and  DIVIDER  seem 
to  be  present  Parameters  were  measured  on  these  beams  (theoretical  beam 
parameters).  Using  the  optimum  detectioa  filter-band  for  GERESS  (lost 
1992a),  i.e.,  1.0  -  3.0  Hz,  a  weak  cmset  for  FLOYDATA  is  visible  (Fig.  2-6). 
Parameters  of  FLOYDATA  were  measured  on  the  corresponding  theoretical 
beam.  For  all  4  events,  the  fk  analyses  produced  less  reliable  results  due  to  the 
low  SNR. 

A  detection  threshold  of  mt,  (NEIC)  =  4.3  is  estimated  using  16  NTS- 
events  recorded  at  GERESS  in  1990-1993.  This  is  done  by  coirqraring  the 
measured  SNR  values  (Table  2-2)  with  the  threshold  of  the  GERESS  STA/LTA 
detector  of  3.8  on  a  test-beam  (e.  g.,  indirect  estimatitm  method;  Ringdal, 
1975).  The  calculated  GERESS  m|,-values  are  in  good  agreement  widi  the  net¬ 
work  magnitudes  published  by  NEIC  GERESS  did  not  detect  (XtSO  with 
m|,  =  4.2  (GERESS).  Also  GALENA  m),  =  4.0  (GERESS)  was  only  detected 
after  reprocessing  at  a  SNR  of  2.1.  On  the  other  hand,  HUNTERS  TROPHY 
with  mi,  =  4.2  (GERESS)  and  also  the  NFT  with  mi,  =  4.0  (GERESS)  were 
clearly  detected  with  a  SNR  above  5.  The  influence  of  noise  (day  /  night;  wind) 
as  w^  as  station  outages  affect  the  SNR  and  can  explain  this  vuialnlity. 

In  order  to  obtain  a  tough  estimate  cm  the  yield  fiom  the  single  site 
GERESS,  we  fit  a  straight  line  tiitough  the  2  known  yields:  The  NPT 
(m^  (GERESS)  =  4.0)  had  an  equivalent  yield  of  1  kt  following  Denny  and 
Zucca  (1993).  The  JVE  (KEARSARGE)  had  a  yield  close  to  ISO  kt  but  pro¬ 
duced  a  body  wave  magnitude  of  only  5.4  (GERESS).  On  the  other  hand, 
explosions  lite  JORNADA  produced  a  magninide  of  5.9  at  a  comparaUe  yield 
(OiqMa  et  aL,  1989).  JORNADA  was  detonated  below  the  water  table  and 
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mformation  from  three-component  and  array  data. 

3.2.3  Determination  of  Polarization  Directions 
For  three-component  data  =  (uiY  with  components  j  =  1,2,3,  the 
covariance  matrix  at  zero-lag  is  computed  as 


where  N  is  the  number  of  samples  in  the  time  window.  Obviously,  S  is 
symmetric,  so  that  its  eigenvalues  are  real  and  nonnegative.  We  choose 
a  set  of  3  orthonormal  eigenvectors  a^  .  A^a^  (no  summation)  then  gives  the 
axis  of  the  polarization  ellipsoid  (see  f.e.  Kanasewich,  1981).  If  the  particle 
motion  is  linearly  polarized,  then  A^  =  A^  =  0.  For  elliptical  particle  motion 
(within  a  polarization  plane),  we  have  A®  =  0.  In  general,  however,  the 
particle  motion  is  ellipsoidal  due  to  presence  of  different  types  of  noise.  K  we 
take  components  (1, 2, 3)  as  (Vertical  up.  North,  East),  the  incidence  angle 
is  =  cos“^lui^|  and  the  backazimuth  is 
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__  (  fan~^(|*^-)  for  <  0 
[x  +  tan“*(||r)  forai^>0 


k  simple  check  of  the  quality  of  the  measurement  is  given  by  the  rectilinearity 


2A» 


(2) 


which  is  equal  to  1  for  linear  polarization  and  0  in  the  worst  case  of  "spheri¬ 
cal”  polarization,  when  all  eigenvalues  have  the  same  size. 

Figme  3.2  shows  a  seismogram  example  demonstrating  a  few  characteris¬ 
tics  of  the  eigenvector  analysis.  For  this  data  example  (year  1992  day  48  hour 
19)  with  an  epicentral  distance  of  385km  and  backazimuth  of  235°  a  number 
of  regional  phases  are  visible.  The  main  phases  for  this  distance  range  are 
Pn,  Pg,  Sn  ond  Sg,  They  are  visible  in  the  three-component  seismogram 
on  the  top.  The  lower  traces  give  polarization  attributes  computed  using 
a  gliding  window  with  a  length  of  2  seconds:  The  root-mean-square  ampli¬ 
tude  of  the  three-component  seismogram  or  "three-component  size”  (fourth 
trace)  X^  =  Sjj  shows  Pg,  Sn  aud  Sg  clearly.  The  first-arriving  Pn 

shows  up  with  somewhat  smaller  three-component  size,  but  with  the  high¬ 
est  rectilinearity  (2),  with  close-to-ideal  linear  polarization.  A  line  gives  the 
critical  rectilinearity  level  of  0.9.  Only  Pn  and  Pg  exceed  this  level.  Since 
the  eigenvalues  are  related  to  energy,  error  bars  can  be  obtained  from  the 
F-test  (see  Bokelmann,  1992).  Sg  shows  up  with  higher  rectilinearity  than 
other  S  phases  and  also  substantial  three-component  size.  The  horizontal- 

vertical  ratio  (H/V)  snd  similarly  the  radial-transverse  ratio  (R/T) 


ddence  angles  within  the  shear  wave  window  and  sources  close  enough  to  the 
receivers  to  not  be  sensitive  to  more  distant  anisotropy.  These  conditions  are 
difficiilt  to  meet  for  our  region-  In  contrast,  for  the  study  of  local  anisotropy 
P-wave  polarization  data  may  be  used  from  local,  regional  and  teleseismic 
events,  if  the  3-component  array  meets  the  measurement  requirements. 
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Figure  3.8:  IllustT&tion  of  plane  P-wave  propagation  in  the  cruat  under  the 
three-component  receiver  array.  In  an  anisotropic  meditun,  the  po¬ 
larization  directions  generally  deviate  from  the  propagation  directions. 
This  deviation,  here  shown  for  the  incidence  angle  69,  is  used  to  infer 
effective  crustal  anisotropy.  Note  that  the  propagation  direction  typi¬ 
cally  deviates  from  the  predicted  propagation  direction,  due  to  source 
mislocation  and  distant  heterogeneity.  With  observations  of  polariza¬ 
tion  and  propagation  directions  distant  effects  can  be  excluded.  Loca' 
disturbances  of  the  wavefront  are  downweighted  by  averaging  over  all 
three-component  stations. 
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onset  if  the  STA/LTA  ratio  for  a  filtered  trial-beam  exceeds  a  given  threshold 
ranging  between  3.S-4.0.  Since  September  1993,  horizontal  beams  have  been 
included  to  inqnove  on  S-type  detections.  The  next  step  of  the  on-line  process¬ 
ing  is  the  transframatitM)  of  a  filtered  data  segment  at  each  onset  time  (derived 
from  the  detection  time)  into  the  frequency-wavenumber  domain.  As  a  result, 
the  slowness  and  back-azimuth  of  dte  pha%  is  determined.  From  the  slowness 
information,  seismic  phases  are  identified.  The  final  step  of  the  on-line  data 
processing  is  the  location  of  events.  The  seismic  phases  as  identified  in  the  fk- 
analysis  are  associated  to  events  in  this  step.  Fiom  the  arrival  time  difference 
of  regional  phases,  the  distance  to  the  epicenter  is  determined  from  the  Jeffreys 
-  Bullen  travel  time  tables  for  regional  seismic  phases.  Togedter  with  a  mean 
back-azimuth,  the  epicenter  locations  of  local  and  regional  events  are  deter¬ 
mined.  Correcting  for  outages,  20  events  have  been  automatically  located  each 
day  on  average  (Table  1-1).  Figure  1-1  shows  the  automatic  event  locations  for 
Dk.  1993  -  May  1994.  Local  magnitivt^g  automatically  determined  at  Bochum 
(e.  g.  Jost,  1993b)  appeared  low  by  at  least  0.5  magnitude  units  compared  to 
values  published  by  other  Euix^iean  observatories.  A  new  amplitude-distance 
curve  has  been  implonented  on  March  1,  1994.  A  detailed  des^ption  will  be 
given  elsewhere. 

Alert  messages  for  strong  teleseismic  and  regional  events  have  been 
automatically  sent  to  NEIC  (Golden,  CO)  and  other  interested  institutions  in 
near  real-time.  44  times,  a  GERESS  alert  message  for  a  teleseismic  onset  has 
been  used  in  the  NEIC  alerts  in  the  reporting  period;  and  19  times  by  the 
EMSC.  For  re^onal  events,  a  fast  earthquake  information  system  (FEIS)  has 
been  develtqwd  (Schweitzer,  Schulte-Theis,  Jost)  and  operated  at  Bochum. 
Starting  from  a  GERESS  aubxnatic  detection  above  a  certain  threshold,  the  rou¬ 
tine  automatically  accesses  the  detectkm  logs  of  the  German  Regional  Seismic 
Netwmk  (GRSN).  Using  all  availatrie  detections  from  diis  network,  the  routine 
proceeds  to  calculate  a  new  location  which  is  subsequently  e-mailed  to 
interested  institutions.  The  time  difference  between  origin  time  and  transmis¬ 
sion  of  an  alert  message  was  42  minutes  for  mrae  than  84  %  of  all  FEIS  alerts 
(Jul  93  -  Oct  93).  The  shortest  time  was  24  minutes,  the  longest  1  hour.  This 
time  delay  was  mainly  determined  by  the  availability  of  communication  lines. 

The  specific  project  of  re-locating  regional  events  by  a  human  analyst  was 
limited  to  1993.  From  January  1994,  analysts  have  routinely  reviewed  regional 
first  arrivals  together  with  teleseismic  events.  For  this  rqwrting  period,  parame¬ 
ter  data  (onset  tinoe,  anqilitude,  period,  azimuth,  and  slowness)  for  S>783  phases 
(i.e.,  54  phases  per  day)  have  been  sent  to  NEIC  (Golden,  (3o),  ISC,  EMSC, 
and  other  institutions. 

GERESS  is  an  open  station  and  resulu  of  the  on-line  processing  have 
been  sent  to  interested  institutions  via  e-mail  (e.g.  Bundesanstalt  fur 
Geowissenschaften  und  Rohstoffe  (BOR),  EMSC,  Frankfurt  University, 
Geoforschungszentrum  Potsdam  (GFZ),  Gtifenberg  Array,  ING,  KTB  array, 
LOG,  Leipzig  University,  Stuttgart  University,  NORSAR,  (Xclahoma  Geol.  Sur¬ 
vey,  ORFi^S,  and  seismological  institutes  in  Austria,  CZ,  Hungary,  Slovenija, 
Switzerland;  YKA).  In  addition,  the  On-line  Processing  Diqrlay  Manager  (Jost, 


Mesa  and  Yucca  Flat  do  not  form  independent  clusters.  This  result  agrees  with 
previous  observadoos  by  Murphy  et  al.  (1989)  who  found  that  the  seismic 
source  characteristics  of  both  testing  areas  are  identical  fOT  tests  below  the 
water  table.  In  general,  we  observe  large  correlation  coefficients  for  events  at 
Pahute  Mesa  (Hg.  2-9).  There  the  depth  of  burial  is  larger  conqrared  to  events 
at  Yucca  Flat  in^tlying  a  largo’  yield.  The  lower  correlation  values  for  most  of 
the  shots  at  Yucca  Flat  are  probably  due  to  their  small  size  and  hence  the  small 
SNR  of  the  beams  entering  the  asrelatitHi  analysis.  Another  reason  may  be  the 
variability  of  the  geologic  structure  at  Yucca  Flat  (e.g.,  McLaughlin  et  al., 
1987;  Cormier,  1987;  Ferguson,  1988). 


2.4  TUAMOTU  ARCHIPELAGO 

In  the  epicentral  distance  range  between  143°  and  147°,  the  structure  of 
the  Earth’s  core  aitqtlifies  -  like  a  focusing  lois  -  seismic  waves  which  travel 
through  the  Earth’s  core  as  PKP-phases.  The  ctxnpariscHi  of  the  daily  detection 
lists  of  GERESS  with  bulletins  of  international  seismic  data  centers  shows  that 
GERESS  is  therefore  very  sensitive  fen:  the  seismicity  in  the  southern  Pacific. 
The  two  French  nuclear  test  sites  Mururoa  and  Fangataufa  are  located  very 
close  to  each  other  in  the  Tuamotu  Archipelago  (Table  2-1).  They  have  an  epi¬ 
central  distance  to  GERESS  of  about  145.5°,  right  where  the  maximum  of  the 
PKP-caustic  is  observed.  Therefore,  we  expect  a  low  detection  level  for  seismic 
events  at  these  test  sites. 

The  wavefonns  of  14  presumed  Freirch  nuclear  tests  in  1988-1991  are 
shown  in  Hgure  2-10  (optimum  beams,  single  traces  for  events  in  1988,  all 
filtered  0.5  -  2.5  Hz).  events  are  odered  with  reqrect  to  increasing  yield 
estimate  (Table  2-4).  Table  2-4  gives  the  coneqxmding  parameters  measured 
<m  the  optimum  beams  (events  in  1990  and  1991)  or  tm  the  single  traces  for  the 
events  in  1988  (recorded  at  station  HAID).  The  Iiench  authorities  have  not  yet 
published  any  parameters  of  their  tests.  Therefore,  the  start  time  of  the 
waveform  for  the  event  ra  Nov.  14,  1990  was  chosen  to  10  seconds  before  the 
measured  onset  time.  The  other  traces  were  aligned  to  its  first  mininaum.  In 
1991,  c»e  of  the  presumed  explosions  was  not  detected  by  GERESS  (Jul  5, 
1991;  18:(X)).  This  event  was  very  small  and  was  only  detected  by  a  T-phase 
observation  at  Raratonga  (RAR).  The  yield  was  estimated  by  DSIR  (Smitii, 
pers.  comnnmicatiOT)  as  less  than  0.3  kt  Hgure  2-10  also  di^lays  the 
GERESS  theoretical  beam  at  the  expected  onset  time  for  this  event,  but  no  con¬ 
clusive  signal  could  be  identified. 

The  body  wave  magnitude  mi,  is  not  defined  for  events  with  an  qncentral 
distance  of  more  titan  about  100°  and  no  direct  determination  is  possible  for 
these  events  from  GERESS.  Schlittenhardt  (1988a)  introduced  a  body  wave 
magnitude  for  PKP-observations  of  the  French  explosioos  at  the  Otifenbetg 
array.  Although  the  distance  between  the  refereoce  station  Al  (rf  GRF  and 
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the  lower  hemisphere.  For  example  polarization  anomalies  from  Northern 
Italy  and  Switzerland  show  opposite  tangential  effects.  However,  based  on 
polarization  data  alone,  this  behaviour  can  not  be  interpreted  in  terms  of 
heterogeneity  or  anisotropy.  Consistently  occuring  source  mislocation  effects 
for  these  regional  distances  might  also  cause  such  variation.  While  this  may 
be  a  problem  of  either  polarization  vector  data  or  propagation  vector  data 
taken  by  itself,  the  difference  is  insensitive  to  source  location  errors  and  nei¬ 
ther  to  distsmt  lateral  heterogeneity.  Instrument  calibrations  were  performed 
daily.  Substantial  deviations  from  the  nominal  response  were  rare  and  the 
correction  had  a  rather  small  effect  on  the  lower  hemisphere  representation. 

3.2.5  Free  Surface 

Conversion  at  the  free  surface  perttirbs  the  incidence  angle  of  the  polar¬ 
ization  direction.  The  apparent  incidence  angle  for  a  P-wave  is  given 
by 


with  the  velocity  ratio  A  =  Vp/v,  (Aki  and  Richards,  1981).  The  dependence 
of  on  $  is  shown  in  figure  3.4  for  a  number  of  Poisson’s  ratios  /c,  which 
are  related  to  the  velocity  ratio  A  by  /c  =  three  values  are 

K  =  0.2  (A  =  1.63;  dotted  line)  ,/c  =  0.25  (A  =  1.73;  solid  line)  and  k  =  0.3 
(A  =  1.87;  broken  line).  Figure  3.4  shows  that  the  free-surface  effect  causes 
incidence  angle  deviations  of  several  degrees.  The  free-surface  correction 
is  particularly  important  for  large  incidence  angles.  Variation  of  k  from 
K  =  0.25  to  0.2  or  0.3  produces  offsets  within  5®  in  the  incidence  range 
of  our  study  (4®  tc  40®),  The  deviation  between  the  curves  is  within  a 
few  degrees.  For  the  free-surface  correction  (figure  3.3)  we  used  a  crustal 
Poisson  ratio  of  k  =  0.25.  This  corresp>onds  closely  to  the  upper  crustal 
mean  value  observed  in  the  KTB  drill  bole,  which  is  located  about  150km 
northwest  (figure  3.1a)  in  a  similar  geological  setting  (figure  3.1).  Without 
this  correction,  the  lower  hemisphere  has  a  distinctly  different  appearance 
(figure  3.5)  with  a  clear  radial  component.  The  mean  of  this  radial  component 
is  5.1®  oriented  outwards.  Clearly,  the  free-surface  correction  is  an  important 
and  observable  effect,  which  must  be  taken  into  account.  On  the  other  hand, 
we  have  seen  that  there  are  clear  effects  in  figure  3.3  beyond  the  simple  effect 
of  crustal  Poisson  ratio  alone.  For  the  Southern  part  of  figure  3.3  the  radisJ 
effect  almost  vanishes. 

3.2.6  Array  Measurements 

Propagation  directions  (^«,4j)  can  be  extracted  from  array  data  either 
from  travel  times  or  from  fk-analysis.  Characteristics  of  this  problem  have 
been  studied  by  Bokelmann  (1993),  where  it  was  shown  that  small-array  data 
like  the  ones  in  this  paper  can  produce  propagation  vector  data  with  reason¬ 
able  imcertainty  given  that  topography  (figure  3.1)  is  taken  into  accotmt.  We 
obtsdn  the  mislocation  pattern  of  figiue  3.6.  Conversion  from  slowness  p  to 
incidence  an^  4.  =  sin~^{pv)  was  done  uring  the  local  velocity  estimate  of 
v=5.2  km/s  obtained  from  3D  array  analysis  (Bokelmann,  1993).  In  that  pa- 
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The  Saturated  Crack  Model  (Lower  Hemisphere  Projections) 


Fieure  3  10-  Lower  hemisphere  representation  of  P-wave  velocity  and  po¬ 
larization  anomalies  for  (r?,  r)  =  (0.82,  1.00)  corresponding  to  the 
saturated-crack  model  HCSl  of  Crampin  (1984). 
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1991b)  has  been  upgraded  several  times  to  make  the  diverse  local  data-bases 
available  for  near  time  information  on  parameter  data  (Fig.  1-2).  The  Data 
Request  Manager  (DRM)  by  Krake  Inc.  (under  contract  fiom  BGR,  Hannover) 
has  been  operated  in  trial  mode.  Its  performance  has  remained  insufficient  and 
a  firequently  announced  upgrade  has  been  waited  for.  Eventually,  the  DRM 
should  enable  the  transmission  of  data  (GSE,  SEED)  and  interface  to  the  Ger¬ 
man  Regional  Netwcnk  of  Broad-Band  Stations. 


Table  1-1 

Statistics  on  the  GERESS  Bochum  On-line  Processing  fOT  Dec  1993-May  1994 
_ number  number/day 


detections 

63608 

349 

f-k  analyses 

50531 

278 

locadons 

3690 

20 

1.3  GERESS  Array  Hub 

The  GERESS  array  has  shown  stable  tqreradon  fiom  Dec.  1993  •  May 
1994.  Table  1-2  shows  die  uptime  of  the  system  including  the  Bochum  data 
center  (uptimes  directly  conespond  m  the  data  archived  in  Bochum).  The  aver¬ 
age  uptime  was  98.9  %,  i.e.,  the  total  outage  amounts  to  1  day  23  hours.  The 
highest  uptime  (99.99  %)  was  observed  in  March  1994;  the  lowest  uptime  (9S.3 
%)  occurred  in  Dec.  1993  due  to  the  sluggish  response  of  die  German  Telekom 
in  replacing  a  data  modem  dannaged  by  thunderstorm.  Table  1-3  further 
specifies  the  causes  of  the  obsoved  outages.  The  technical  status  of  GERESS 
is  summarized  in  chronological  order  in  Appendix  1-1.  Results  of  a  recent 
instrument  cafitvation  are  given  in  Appendix  1-2. 

Tabic  1-2 

GERESS  uptime  Dec.  1993  -  May  1994  (incl.  Bochum  data  center) 


Month  uptime 
% 


Dec 

95.30 

Jan 

99.95 

Feb 

99.96 

Mar 

99.99 

Apr 

99.96 

May 

98.42 

17 


GERESS  is  only  about  200  km.  we  cannot  use  his  results  to  estimate  m|,-values 
because  the  site  conditions  differ  substantially.  The  an^tlitude-^tance 
behavior  of  the  PKP-phases  varies  drastically  in  the  distance  range  near  145.5° 
(e.  g.,  Houard  et  al.,  1993).  and  we  corrected  the  observed  log  (.ATT)  values 
with  the  attenuatitm  curve  of  Blandfotd  and  Sweetser  (1973)  fm  PKP-phases. 
\^th  the  correction  value  of  3.50  fiff  an  epicentral  distance  of  145.5°.  we  deter¬ 
mined  the  nib(PKP)-values  for  the  presumed  French  tests.  For  11  of  these 
events,  our  rttt,(PKP)-values  can  be  compared  with  m|,-determinations  by 
NEIC.  Although  our  m{,  values  are  slightly  smaller  than  the  network  magni¬ 
tudes  by  NEIC.  the  general  agreement  betwemi  the  m|,-values  calculated  with 
different  methods  confirms  the  proposed  usage  of  PKP-phases  for  regular 
world-wide  monitoring  (Harjes,  1985). 

Figure  2-11  shows  the  estimated  m),(PKP)-values  for  GERESS  plotted 
against  the  yields  estimated  by  DSIR  (Table  2-5,  Smith,  pers.  communication). 
The  line  follows  the  linear  relationship  Qeast  squares  fit)  between  log(yield)  and 
mj,  (PKP)-values: 

logio  (yield)  =  1.24  mb  (PKP)  -  4.88  .  (2) 

This  relation  is  very  similar  to  the  results  of  Schlittenhardt  (1988a)  who  got  for 
mb(PKP)  based  on  GRF  measurements; 

logio  (yield)  =  1.12  itH,  (PKP)ow:  -  4.58  .  (3) 

The  two  solid  lines  in  Hgure  2-11  give  the  iqrper  limits  for  die  smallest  not 
observed  explosion.  On  average,  the  GERESS  detection  threshold  is  about 
mb(PKP)  =  3.8  cone^Kniding  to  a  yield  estimate  of  0.7  kt  For  the  network 
(NEIC)  magnitude,  the  GERESS  detectim  threshold  is  4.0. 

Nuclear  tests  in  die  Tuamotu  Archipelago  have  been  rqiorted  to  be 
detonated  either  on  the  island  of  Muniroa  or  on  the  island  of  Fangataufa  with 
about  40  km  in  between.  The  cotiespcmding  variadcHi  in  slowness  and  azimuth 
for  PKP-onsets  is  so  small  that  even  large  aperture  arrays  (e.g.,  GRF)  are  not 
able  to  discriminate  between  the  2  test  sites.  Therefore,  we  used  a  different 
approach  and  performed  a  cross-conelatkm  of  all  recorded  events  (1988-1991). 
We  used  a  30  second  time  window  of  the  tqitimum  GERESS  beams  filtered 
between  0.5-2.5  Hz.  Figure  2-12  shows  a  cross-conelation  sequence  of  die 
event  on  Nov.  21,  1990  with  all  events  in  1990-1991.  Note  t^  the  cross- 
correlation  coefficients  of  the  three  evmts  on  the  bottom  of  the  figure,  rqiott- 
edly  placed  (»  Fangataufa,  are  significandy  different  (very  small)  fiom  those 
reportedly  placed  on  Muniroa.  The  cross-correlation  coefficients  are  also  small 
for  the  fim  two  events.  The  first  evmt  has  been  classified  as  not  observed,  and 
the  second  (yield  in  the  order  of  1  kt),  although  observed,  is  contaminated  with 
too  much  noise.  On  the  other  hand,  the  cross-correlation  sequence  indicates 
that  the  observed  waveform  is  from  the  cotieqionding  event  at  Tuamotu.  Fig¬ 
ure  2-13  shows  the  cross-conelatioa  matrix  obtained  by  using  the  maximum 
cross-conelation  value  in  the  filter-band  0.5-2.5  Hz  and  the  event  numbers  from 
Table  2-5.  Table  2-6  lists  the  minimum  and  muTimum  correlation  coefficients 
for  each  event  with  respect  to  the  Muniroa  and  Fangataufa  clusters.  Rqieating 
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Figure  2-11:  Relaticm  of 
Tuamotu  Archipelago. 
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per,  the  reader  also  finds  a  discussion  on  the  consistency  of  propagation  vector 
estimates  from  examples  of  Nevada  and  IViamotu  nuclear  events.  These  data 
are  included  in  this  data  set.  In  figure  3.6  they  show  comparatively  small 
deviation  documenting  consistency  for  the  array  data.  Nevertheless  these 
deviations  showed  up  as  statistically  significant  (Bokelmann,  1993).  There 
are  even  larger  anomalies  in  these  data  thotigh.  The  Northern  deviation  for 
steeply  incident  phases  from  the  North  which  was  observed  for  polarization 
data  (figure  3.3)  is  not  apparent  in  the  array  data.  Interestingly  the  pattern 
for  Yugoslavia  and  Romania  are  very  similar  to  observed  polarization  anoma¬ 
lies  in  figure  3.3  as  well  as  the  rapid  variation  of  Switzerland  and  Northern 
Italy.  This  suggests  that  either  lateral  heterogeneity  or  source  mislocation  is 
affecting  propagation  and  polarization  directions. 

For  an  epicentral  distance  of  about  4001cm  the  mean  tangential  deviation 
between  Switzerland  and  Northern  Italian  events  would  require  a  relative 
source  mislocation  of  about  70km  occuring  consistently  in  the  PDF  solutions 
(Preliminary  Determination  of  Epicenters)  between  the  two  source  regions. 
While  this  is  not  impossible,  lateral  heterogeneity,  the  Alpine  root,  might  also 
cause  the  2izimuthal  deviation,  or  respectively,  the  relative  source  mislocation. 
Both  would  simultaneously  perturb  propagation  and  polarization  directions. 

Nevertheless,  the  procedure  of  relative  analysis  of  polarization  and  prop>- 
agation  data  gives  data  which  are  independent  of  the  source  location  and 
distant  lateral  heterogeneity.  The  relative  data  (polarization  -  propagation 
vectors)  from  the  94  events  are  shown  in  figure  3.7.  However,  Nevada,  Tu- 
amotu,  and  other  regions  still  show  up  with  consistent  effects.  Naturally,  the 
random  scatter  is  larger  for  these  differences  than  for  the  individual  data.  For 
two  events  the  difference  data  were  clearly  dominated  by  noise  and  had  to 
be  removed.  It  may  seem  surprising  that  source  mislocation  and  lateral  dis¬ 
tant  heterogeneity  don’t  affect  the  relative  data.  In  fact,  the  method  doesn’t 
require  knowledge  of  source  locations  at  all.  The  propagation  vector  is  used 
to  determine  the  location  in  the  lower  hemisphere.  Measured  from  this  di¬ 
rection,  the  polarization  vector  deviation  indicated  presence  of  anisotropy. 
This  allows  studying  anisotropy  under  the  receivers.  Lateron,  we  use  the 
data  of  figure  3.7  in  an  inversion  for  anisotropy.  But  first,  we  review  wave 
propagation  theory  in  anisotropic  media. 

3.3  REVIEW  OF  PROPAGATION  IN  ANISOTROPIC  MEDIA 


Here  we  wish  to  merely  state  a  few  features  of  partictilar  importance  for 
this  study.  Comprehensive  treatment  of  the  topic  has  been  given  f.e.  by 
Crampin  (1981). 

In  a  general  anisotropic  medium,  the  solution  for  a  propagating  plane  wave 
front  follows  from  the  Ansatz 

u(r,  t)  =  a(r)/(t  -  nr/ c)  (4) 

where  a(r)  describes  amplitude  and  polarization  of  the  wave,  while  the  time- 
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Table  1-3 

Causes  of  GERESS  downtimes  Dec.  1993  -  May  1994  (incl.  Bochum  data  center) 
(The  total  downtime  was  171922  sec  «  id  23h  45m  22  sec) 


cause 

downtiine 

downtime 

[sec] 

[%1 

Bochum  data  acq.  software 

17075 

9.9 

Bochum  data  acq.  software  (maintenance) 

0 

0 

Bochum  data  acq.  workstation 

5 

0 

Bochum  CIM 

0 

0 

Bochum  power  outage 

41726 

24.3 

total  outage  Bochum 

58806 

34.2 

Communication  line 

75685 

HUB  lAC  crashes  (1) 

98 

0.1 

HUB  thunderstonn  damage 

19774 

11.5 

HUB  maintenance 

3659 

2.1 

HUB  power  outage 

13769 

8.0 

total  outage  HUB 

37300 

21.7 

other  causes 

131 

0.1 

GERESS  has  used  data  acquisition  equipment  from  Teledyne-Geotech  (i.e. 
RDAS-200,  lAC)  which  is  no  longer  in  production.  GERESS  is  a  prototype 
installation  and  inherent  problems  have  remained  (e.g..  Golden  et  1991). 
Due  to  the  severity  of  th^  problems,  we  will  continue  to  list  them  to  remind 
the  community  of  possible  problems  inherent  in  GERESS  data. 

-  The  RDAS-2(X)  occasionally  sends  electronic  noise  instead  of  seismic  data 

after  power  iq)  or  self-reset  This  state  can  be  cleared  only  by  manual 
intervention  (reset). 

-  If  the  Ipps  timing-signal  is  bad  (optical  modems,  fiber  <q>tic  link),  or  the 

clock-board  in  the  RDAS-200  deteriorates  before  completely  failing,  or 
after  any  self-reset  or  after  any  power  cycle,  data  can  te  desynchronized 
by  one  second  or  by  less  than  a  second  (mostly  0.2  seconds).  Only  an 
array  installation  of  RDAS-2(X)’s  has  the  chance  to  detect  this.  Each  day,  a 
careful  inspection  GERESS  records  has  been  performed  to  insure  pnqter 
synchronizjttion.  After  a  desynchronization  hu  been  detected,  one  or 
several  manual  power  cycks  usually  synchronizes  the  station. 


-  The  RDAS-2(X)  is  very  sensitive  to  spikes  in  input  voltage  (thunderstonns) 
causing  outages  or  s^-resets. 
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the  conelaticMi  analysis  including  all  lecoided  events  in  1988  shows  that  the 
difference  of  cocrelation  coefficients  between  those  two  clusten  decreases.  We 
e^lain  this  by  the  larger  noise  level  of  the  single  station  recoids  (evoits  in 
1988).  Other  pass-bands  (0.7-2.0  Hz,  1.0-3.0  Hz)  genendly  support  tiM  presmt 
finding,  although  the  observed  grouping  of  coneladon-coefificients  is  C9>timum 
for  most  events  in  die  pass-band  of  O.S-2.5  Hz.  The  cross-correlation  sequence 
of  the  event  cm  May  25,  1988  (#  11)  which  had  the  largest  repented  yield  on 
Muruioa  (83  kt)  also  agrees  with  the  present  finding,  where  an  imexplained  de- 
ccnreladcm  with  die  event  cm  June  2,  1990  (#  8)  is  observed  (Hg.  2-13).  Hgure 
2-14  shows  the  dendrogram  based  cm  a  sin^  linkage  clustering  analysis.  We 
clearly  see  the  different  Fangataufa  and  Mururoa  clusters.  The  event  #  2 
(07/05/1991;  1  kt)  has  a  very  small  SNR  which  makes  a  reascmable  analysis 
inqmssible.  On  the  other  hand,  we  see  some  varial^ty  in  the  Mururoa  cluster, 
e.g.  events  #  4  (18/5/1991),  #  10  (21/11/1990),  and  #  9  (15/07/1991)  ^ipear  to 
form  a  subgroup  and  also  the  events  #  6  (11/05/1988),  #  5  (04/07/1990),  and  # 
7  (14/06/1991).  This  observation  may  result  firom  the  observed  local  site 
effects  (Crusem  and  Caristan,  1992). 

Following  Guille  et  al.  (1993),  the  atoll  of  Mururoa  (28  km  by  10  km;  age 
11.8-10.6  My)  is  part  of  the  volcanic  chain  feom  the  hot-spot  near  the  island  of 
Pitcairn.  Mururoa  sits  on  ic^  of  a  submarine  plateau  of  130  km  by  30  km 
oriented  N80E.  The  existence  of  this  plateau  is  intimately  linked  to  the  Austral 
feacture  zone  which  explaitts  die  atypical  orientation  of  this  plateau.  Although 
the  island  of  Fangataufa  (10  km  by  5  kzn;  age  11.5-9.6  My)  is  situated  only  40 
km  soudi-east  of  Mururoa,  the  influence  of  this  firacture  Zone  is  minor.  The 
morphology  of  die  basonent  is  star  studied,  quite  differmt  from  the  elcmgated 
plateau  of  Mururoa.  Typical  tectonic  directions  on  Fangataufa  follow  those  of 
the  southern  Pacific  (NS,  N80E,  and  N115E).  The  bathymetry  also  supports  a 
generic  difference  between  both  islands:  the  strait  separating  bodi  islands  shows 
water-depths  larger  than  3  km  for  about  20  km.  As  a  result,  both  islands  exhi¬ 
bit  different  tectonic  settings  that  can  be  re^xmsible  for  the  observed  clustering. 

We  ctmclude  that  based  on  the  limircrf  data  set  (14  events),  the  optimum 
GERESS  beams  can  be  used  in  a  correlaticm  analysis  (waveform  matching)  to 
successfully  sqiarate  the  Mururoa  and  Fangataufa  clusters.  French  authorities 
have  not  yet  published  any  infotmatitm  on  their  tests.  Therefore  we  cannot 
know  to  which  degree  the  above  results  r^ect  different  device  types. 
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dependence  is  given  by  /()  depending  on  the  phase  velocity  c.  FVom  the 
equation  of  motion,  we  get 

muai  =  c^ai  (5) 

with  the  Christoffel  matrix 

mu  =  CijktnjTtk/p  (6) 

(6)  has  the  form  of  an  eigenvalue  equation  with  the  square  of  the  phase  ve¬ 
locity  appearing  as  eigenvalue.  Hence  polarization  vectors  a  serve  as  eigen¬ 
vectors.  Due  to  the  symmetry  of  the  elasticity  tensor  Cjjtj,  we  can  always 
choose  a  set  of  3  mutually  perpendicular  eigenvectors.  Below  we  will  give 
examples  of 

The  typical  way  to  solve  (6)  is  to  specify  a  propagation  direction  n  and 
solve  for  eigenvalue  and  eigenvector  for  each  of  the  three  polarizations  (in 
the  nondegenerate  case).  It  is  clear  that  in  general  propagation  and  polar¬ 
ization  vectors  don’t  coincide  for  waves  propagating  in  anisotropic  media. 
For  isotropic  media  that  was  the  case.  Therefore  we  can  use  such  deviations 
to  determine  effective  medixun  anisotropy. 

For  the  case  of  P-waves,  we  will  discuss  this  discrepancy  more  closely 
below,  since  both  of  these  directions  can  be  observed  using  arrays  and  three- 
component  stations. 

The  wavefront  propagates  in  the  normal  direction  with  the  phase  velocity 
c  =  Vt  is  the  propagation  vector,  which  can  be  estimated  from  data 

of  small  arrays  (see  also  Helbig,  1958). 

The  polarization  vector  a  gives  the  direction  of  local  particle  motion.  It 
is  measured  directly  by  polarization  analysis.  Neither  of  these  two  vector 
should  be  confused  with  the  group  velocity  vector,  which  is  the  direction, 
in  which  energy  propagates  (Musgrave,  1970).  Crampin,  Stephen  and  Mc- 
Gonigle  (1982)  showed  that  for  weak  anisotropy  the  difference  between  po¬ 
larization  and  group  velocity  directions  are  small.  This  is  intuitively  clear, 
since  the  energy  should  propagate  in  a  direction  near  the  particle  motion. 

While  the  angular  deviation  between  phase  and  group  velocity  vectors 
received  quite  some  interest  in  the  literature  (Crampin,  1981;  Crampin  et 
al.,  1982),  relation  between  phase  velodfy  and  polarization  vectors  has  been 
discussed  for  model  crystals  (Crampin,  1978,  1981;  Crampin  et  al.,  1982), 
but  due  to  a  lack  of  appropriate  data,  there  are  only  few  applications  to  real 
data  so  far  (White  et  al.,  1983;  Li,  Leary  and  Aid,  1987;  de  Parscau,  1991a). 
However,  deviations  of  polarization  directions  frt>m  propagation  directions 
can  be  substantial  as  we  will  see  below.  In  the  following  we  want  to  illustrate 
this  effect  and  lateron  exploit  it  in  an  inversion  for  effectively  anisotropic 
structure.  It  is  important  to  note  that  the  polarization  depends  only  on 
local  anisotropy.  This  is  illustrated  in  figure  3.8. 

In  the  context  of  anisotropy,  P-wave  polarizations  have  been  studied  for 
VSP  borehole  data  by  de  Parscau  (1991a)  and  Li,  Leary  and  Aid  (1987). 
However,  with  a  ID  array  in  the  borehole,  there  is  no  control  over  off-azimuth 
arrivals. 
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TABLE  3.1:  Events  Used  in  this  Study 

(Location  and  time  from  PDE  Monthly  Btilletin,  if  not  stated  otherwise) 


Year,  Day,  Hour,  Min,  Sec.  Latitude  Longitude  Depth  Phase 


1990  123  01  03  38 

43.315 

19.890 

5. 

Pn 

1990  240  20  21  22 

36.267 

27.218 

39.9 

Pn 

1990  246  10  48  33 

45.915 

15.873 

20. 

Pn 

1990  253  12  17  36 

70.770 

-13.803 

10. 

P 

1990  258  23  07  43 

64.655 

-17.617 

10. 

P 

1990  285  17  30  00 

37.250 

-116.490 

0. 

P 

1990  297  15  04  14 

73.360 

54.670 

0. 

P 

1990  318  18  11  58 

-22.258 

-138.805 

0. 

PKP 

1990  325  16  59  58 

-21.9 

-138.980 

0. 

PKP 

1990  331  04  37  58 

43.853 

16.633 

24. 

Pn 

1990  331  04  51  36 

43.895 

16.641 

10. 

Pn 

1990  347  00  24  26 

37.3 

15.438 

11.1 

Pn 

1990  348  03  21  27 

39.347 

15.355 

276. 

Pn 

1990  350  15  45  41 

41.361 

43.715 

33. 

P 

1990  355  13  12  53 

-18.891 

-177.971 

457. 

PKP 

1991  038  07  12  48 

47.590 

15.503 

10. 

Pn 

1991  042  15  43  44 

44.871 

6.704 

14. 

Pn 

1991  059  15  29  42 

51.690 

16.360 

0. 

Pn 

1991  074  03  24  09 

34.343 

26.389 

7.1 

Pn 

1991  078  02  51  26 

39.260 

20.429 

10. 

Pn 

1991  094  19  00  00 

37.296 

-116.313 

0. 

P,PcP 

1991  100  01  08  40 

37.359 

36.221 

10. 

P 

1991  111  08  51  56 

37.972 

19.952 

39.8 

Pn 

1991  116  22  35  11 

43.446 

16.244 

12. 

Pn 

1991  119  09  12  48 

42.453 

43.673 

17.2 

P 

1991  119  18  30  42 

42.503 

43.899 

14.3 

P 

1991  120  03  40  36 

51.690 

16.206 

10. 

Pn 

1991  122  10  15  17 

47.929 

16.209 

10. 

Pn 

1991  123  23  41  02 

42.647 

43.263 

10.8 

P 

1991  136  02  06  17 

52.309 

7.649 

52.4 

Pn 

1991  138  17  14  59 

-21.832 

-139.014 

0. 

PKP 

1991  143  19  42  56 

51.416 

15.850 

10. 

Pn 

1991  149  18  59  58 

-22.256 

-138.794 

0. 

PKP 

1991  149  20  24  40 

45.016 

8.213 

10. 

Pn 

1991  154  10  22  40 

40.048 

42.859 

27.8 

P 

1991  165  17  59  58 

-21.944 

-138.988 

0. 

PKP 

1991  166  00  59  20 

42.461 

44.009 

9.4 

P 

1991  185  06  26  32 

42.387 

44.116 

20.4 

P 

1991  191  23  57  20 

51.424 

16.217 

0. 

Pn 

1991  196  18  09  58 

-21.877 

-138.963 

0. 

PKP 

1991  200  01  19  52 

45.344 

21.123 

10. 

Pn 

88 
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Figure  3.12:  Fitting  of  the  saturated-crack  model:  a)  misfit  in  degree^ 
depending  on  symetry  axis  orientation  from  O'to  90°incidence,  b) 
best-fitting  model  shown  in  a  partial  lower  hemisphere  display  (0®to 
55°incidence).  The  symmetry  axis  has  strike  208°  and  dip  36°  to  the 
South.  The  lowest  contour  line  gives  the  critical  level.  Polarization 
deviations  are  given  for  0°,  10°,...,  50°incidence. 
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-  The  RDAS-200  occasionally  soids  data  with  corrupted  tiine  stamps  in  the 

future. 

-  Meteorological  data  are  not  availaUe  for  GERESS. 

-  GERESS  runs  on  DCF-77. 

At  Bochum  data  center,  we  have  continued  to  keep  GERESS  operational  to 
jMovide  the  community  with  useful  data.  We  hope  to  gather  sufficient  q>are 
parts  to  keep  the  current  hardware  running  as  long  as  practical.  For  goreral 
scientific  purposes,  this  situation  appears  acceptaUe.  But  for  the  special 
requirements  in  a  treaty  situation,  we  recommend  an  upgrade  of  the  present 
hardware  to  a  nxxe  robust  system. 

1.4  Conclusion 

GERESS  is  the  most  sensitive  seismological  station  in  Central  Europe  for 
regional  events.  The  teleseismic  detection  threshold  of  GERESS  ranges  glo¬ 
bally  between  m|,  4.0  -  4.S  (this  report). 

At  the  GERESS  data  center  Bochum,  the  observatory  routine  has  included 
dau  archival  and  automatic  on-line  analysis.  From  Dec.  1993  to  May  1994.  20 
regional  events  have  been  automatically  located  each  day  on  average.  The 
GERESS  data  center  Bochum  is  die  only  German  stttion  that  provides  NEIC 
and  interested  institutions  with  alert  messages  (teleseismic  and  regional). 

In  Bochum,  extensive  parameter  data  for  54  regional  and  teleseismic 
phases  have  been  interactively  reviewed  each  day  on  average.  Results  of  this 
re-analysis  GERESS  data  have  been  used  NEIC  ISC  EMSC,  and  other 
interested  institutions.  The  Bochum  data  center  is  an  open  station  similar  to  the 
German  Regional  Network  of  Broad-Band  Stations  (GRSN).  In  addition, 
Bochum  serves  as  data  center  where  various  bulletins  (e.g.,  alerts  from  NEIC, 
EMSC  NORSAR,  SED,  and  ING;  onset  lists  from  GRSN,  Austrian  stations, 
Czech  stations,  and  YKA;  locations  from  NEIC  (C^D,  FDE  (weekly),  PDE 
(monthly)),  IDC  (GSETT-3),  K3N,  IMS,  LDG,  and  SED;  fast  CMT’s  from  Har¬ 
vard,  USGS,  and  Tokyo  Univenity)  are  available  for  the  interested  user. 

The  GERESS  dau  center  Bochum  has  aided  the  CSS  to  lest  the  experi¬ 
mental  GSETT-3  International  Data  Crater  (IDC  VO,  since  May  1993).  Within 
this  concqit,  GERESS  served  as  a-station  and  acted  as  a  node  to  access  data 
from  the  GRSN,  serving  as  ^-stations. 

The  GERESS  data  center  Bochum  has  been  leqxmsiUe  for  GERESS 
maintenance.  The  average  GERESS  uptime  which  includes  the  Bochum  data 
center  has  increased  to  99.65  %  during  the  first  5  months  of  1994.  For  com¬ 
parison,  GERESS  uptimes  were  912  %  in  1993,  96.9  %  in  1992,  and  91.5  % 
in  1991.  This  again  demonstrates  the  high  level  of  reliability  of  the  installation, 
its  successful  maintenance,  and  the  high  motivation  of  key  personnel 
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2.5  LOP  NOR 

GERESS  recoided  S  of  the  6  presumed  nuclear  explosions  at  die  Chinese 
test  site  Lop  Nor  in  1990-1994.  Figure  2-lS  displays  the  optimum  beams 
filtered  0.5  -  2.5  Hz.  Parameters  measured  cm  the  c^dmum  beams  are  given  in 
Table  2-7.  The  observed  body  wave  magnitudes  are  unexpected  low  compared 
to  the  network  magninides  (NEIC).  This  observed  discrepancy  suggests  that 
GERESS  has  a  signifirant  station  residual  for  Nop  Lor  events.  However, 
GERESS  generally  shows  a  negative  magnitude  residual  for  events  in  western 
China  and  the  Ifimalaya  region  of  -0.51  +/-  0.26  magnitude  units.  With  a 
STA/LTA  detection  threshold  of  3.8  of  the  GERESS  detector,  this  implies  a 
detectitHi  threshold  for  events  at  Lop  Not  of  nO),  (NEIC)  =  4.6. 

For  4  events,  surface  wave  magnitudes  could  be  measured  at  GERESS 
(Table  2-7).  However  for  the  events  in  1990  and  1993,  the  magnitude  values 
give  an  upper  bound  only  since  the  traces  were  disturbed  by  interfering  noise. 
The  classical  ooi,  /  MS-ratios  at  GERESS  cleariy  ctmfirm  the  explosion  source 
type  for  all  4  events.  On  the  othOT  hand,  Lc^  Nor  records  show  a  compara¬ 
tively  strong  Love-wave  radiatiOTi  towards  European  staticms,  which  is  not 
expeoed  for  an  explositm  type  source.  Note  that  die  surface  wave  magnitude 
measured  on  the  horizontal  conqxxient  instruments  is  significantly  larger  than 
the  magnitiiA;  measured  on  the  vertical  instrument  (Table  2-7)  because  the 
Love  wave  shows  the  largest  amplitude  in  the  surface  wave  group. 

The  test  on  May  21,  1992  was  tiie  largest  underground  explosion  sriiich 
has  ever  been  observed  with  broad-band  statitnis  world-wide.  To  demonstrate 
the  sinqtle  pulse  foim  of  relatively  high  frequency  for  this  large  event,  Hgure 
2-16  shows  the  deccmvolved  P  and  PcP  (xisets.  In  both  cases,  a  long  periodic 
down  swing  after  the  first  onset  is  seen  which  can  be  addressed  as  an  effect  of 
tectonic  release  related  witii  this  event  (Gao  and  Wallace,  1992). 


2.6  NOVAYAZEMLYA 

On  October  24,  1990,  the  Soviet  Union  carried  out  its  last  nuclear  test 
(Table  2-8).  This  test  was  cOTiducted  on  Novaya  Zemlya  at  a  distance  of  30.4° 
to  GERESS.  The  records  are  dominated  by  the  very  prominent  P-wave  coda 
(Hg.  2-17,  optimum  beam,  filtered  0.5  -  2.5  Hz),  and  a  very  conqtlex  pulse  of 
the  direct  P-wave.  This  feature  must  be  etqtlained  by  contributions  fiom  die 
direct  waves  the  upper  mantle  triplications,  by  converted  phases,  multipath- 
ing,  and  by  a  lateral  heterogeneous  structure  of  crust  and  upper  mantle  in 
Europe.  The  obsenred  magnittide  ni|,  «  5.5  for  this  explosion  is  in  agreement 
widi  NEIC  (tn^  a  5.7).  The  classical  mi,  /  MS-ratio  (5.5  /  4.4  =  1.3)  at  GERESS 
cleariy  confirms  the  explosion  source,  ^tfa  the  observed  SNR  of  90.3  the 
detection  threshold  at  GERESS  can  be  estimated  at  m|,  (NEIC)  «  4.3. 
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Figure  2-13:  CcKielaticm  matrix  of  the  events  firm  Tuamotu  Archipelago  (fcx" 
event  numbers,  see  Table  2-5). 
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We  have  seen  that  our  method  of  relative  analysis  is  insensitive  against 
distant  heterogeneity.  But  how  about  near-receiver  heterogeneity?  Inspec¬ 
tion  of  polarization  hemispheres  for  individual  three-component  instnunents 
showed  some  station-dependent  scatter  which  is  partly  due  to  deviations 
of  the  instnunent  response  from  nominal  values.  Small-scale  heterogeneity 
near  the  receivers  may  also  cause  waveform  distortions.  Hence,  polariza¬ 
tion  data  from  single  three-component  stations  have  to  be  regarded  with 
some  caution.  Such  data  carry  information  about  the  local  conditions  under 
the  station.  For  our  purpose  we  average  over  all  three-component  stations 
of  the  array,  thereby  obtaining  avers^e  polarization  attributes,  where  array 
stations  are  well-distributed  over  an  aperture  of  4  km.  The  effect  of  near- 
receiver  small-scale  heterogeneity  is  expected  to  be  nearly  uncorrelated  from 
station  to  station  and  we  expect  the  averaging  procedure  to  be  effective  in 
averaging  out  most  of  its  effect  on  the  difference  pattern.  The  effect  from 
larger-scale  heterogeneity  on  the  other  hand,  is  likely  correlated  from  station 
to  station.  We  obtain  a  net  effect  in  the  polarization  vectors,  but  simi¬ 
larly  in  the  propagation  vectors.  The  difference  should  therefore  again  be 
approximately  vmaffected.  This  argument  suggests  that  although  scatter  is 
invariably  present  in  difference  patterns  like  figure  3.7,  the  scatter  is  domi¬ 
nated  by  random  fluctuation  and  not  by  biases  from  discrete  local  scatterers. 
A  very  important  feature  of  the  GERESS  array  is  the  absence  of  a  seismo- 
logically  significant  near-surface  low- velocity  layer,  what  can  be  shown  from 
inspection  of  the  train  of  P-S  conversions. 


3.4  UPPER-CRUSTAL  ANISOTROPY 


Anisotropy  in  the  vicinity  of  the  receivers,  however,  does  affect  the  relative 
data  in  figure  3.7.  Which  depth  range  is  it  sensitive  to?  We  have  seen  that 
polarization  is  a  local  property.  For  seismological  phases  of  finite  wavelength 
this  means  that  the  depth  range  equivalent  to  about  a  wavelength  should 
affect  the  observed  polarization.  For  the  typical  dominant  frequencies  in  this 
study  we  have  wavelengths  of  about  6  km.  Hence,  we  expect  to  estimate 
a  spatial  average  of  anisotropic  properties  of  the  top  few  kilometers  in  the 
upper  crust.  Since  we  average  over  all  three-component  stations  in  the  array, 
we  assiune  that  discrete  heterogendty  within  the  array  can  be  ignored  and 
allow  for  an  anisotropic  but  homogeneous  model. 

In  this  context,  there  are  several  phenomena  which  may  cause  effective 
anisotropy.  Intrinsic  anisotropy  of  minerals,  as  well  as  small-scale  structure 
possessing  certain  symetry  lead  to  directional  dependence  of  macroscopic 
elastic  parameters  and  hence  propagation  velocity.  Sedimentary  layering, 
with  a  potentially  prominent  effect,  is  not  present  here.  However,  cracks, 
firactures  and  rock  foliation  (Babuska  and  Cara,  1991),  if  certain  orientations 
prevail,  can  give  rise  to  substantial  effects.  Since  our  data  are  sensitive  to 
a  spatial  average  over  several  kilometers,  we  expect  macroscopic  anisotropic 
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Using  Nuttli’s  (1986,  1988)  linear  relaticnt  between  nii,  ^g)  and  yield  and 
taking  pnq>er  care  of  the  magnitude  bias  between  P-wave  and  Lg-wave  magni¬ 
tudes,  we  obtained  a  yield  estimate  of  50  kt  firom  die  GERESS  P-wave  magni¬ 
tude.  Schlittenhaidt  (1988b-1992)  determined  the  yield  for  the  Novaya  Zonlya 
test  as  70  kt  based  on  GRF  records  (mj,(GRF)=5.9).  For  GERESS,  we 
estimated  a  detection  threshold  in  the  order  of  3  kt  Note  that  this  value  is 
based  on  (me  single  shot  only. 


2.7  EASTERN  KAZAKHSTAN 


During  the  site  survey  for  GERESS,  6  events  fixim  the  framer  Soviet  test 
site  near  the  city  of  Semipalatinsk  were  recraded.  All  tests  were  conducted  at 
the  Shagan  River  (Bali^ian)  subsection  of  the  testing  area  (e.g.,  Thurber  et  al., 
1993).  We  recorded  the  event  on  May  4,  1988  at  the  tenqxmiy  station  HAID, 
and  the  tests  on  Jun  14  and  Sep  14  (JVE),  1988  at  die  temporary  statirai  GRLB 
(Haijes,  1989).  The  three  events  in  Dec.  1988,  Jan.  and  Feb.  1989  were 
recorded  with  the  tenqxxary  9-element  GERESS  test  array  (Harjes  1990a). 
Figure  2-18  shows  the  sin^e  traces  and  the  optimum  beuns,  reflectively, 
filtered  between  0.5  Hz  -  2.5  Hz.  Excqit  for  the  event  on  June  14,  1988,  mag¬ 
nitudes  measured  at  die  GERESS  ate  cotielatB  well  widi  die  ISC  network  mag¬ 
nitudes  (Table  2-9).  Based  on  these  6  events,  an  average  detection  threshold  of 
mi,  (NEIC)  =  4.7  has  been  determined.  This  estimate  gives  an  iqiper  bound 
only,  since  the  25-element  GERESS  array  shows  a  significandy  better  SNR 
inqirovement  by  beamfotming  than  a  single  station  or  the  9-element  test  array. 
Hence,  the  detection  threshold  for  the  25-element  GERESS  array  can  be 
estimated  to  be  lower  than  4.5  fra  eastern  Kazakhstan.  This  estimate  is  conser¬ 
vative  because  the  comparison  of  the  technical  installation  of  GERESS  and  the 
tempraary  test  array  suggests  a  sensitivity  surpassing  the  expected  Vn  SNR 
infRDvement  Ringdal  (1990)  determined  the  50%  detection  dneshold  for 
NORESS  as  mi,  =  2.7  -  3.3.  The  discrepancy  to  die  estimate  fra  GERESS  is  in 
part  due  to  the  well  documented  focusing  effects  in  die  upper  mande  below  the 
test  site  towards  Scandinavia.  Barker  and  Mur^diy  (19^)  determined  station 
residuals  in  Europe  vdiich  amount  to  more  than  1  magninide  unit 

In  Table  2-10,  we  cotrqiare  various  yield  determinations  fra  the  6 
presumed  tests.  Using  Nuttli’s  (1986,  1987)  linear  relation  between  m|,  (Lg) 
and  yield  and  taking  proper  care  of  the  magnitude  bias  between  P-wave  and 
Lg-wave  magnitnHM,  we  obtained  yield  estimates  firom  GERESS  P-wave  mag- 
nitndes.  Included  are  also  yield  estimates  by  Sdilittenhardt  (1988b-I992)  foom 
GRF  records.  Based  on  die  Lg  wave  at  NORSAR,  Hansen  et  aL  (1990)  calcu¬ 
lated  a  yield  of  110  kt  fra  the  JVE  on  Sq>.  14,  1988.  This  value  agrees  with 
oar  estimate  of  108  kt  at  CXRESS.  Rom  hydrodynamic  measurements,  the 
yield  was  determined  from  US  and  USSR  scientists  at  115  kt  and  122  kt, 
reflectively  (Sykes  and  Ekstrom,  1989).  Ringdal  et  aL  (1992)  and  also 


Table  2*10:  Magnitude  Yield  Relations  for  Kazakhstan  from 
6  Events  Recorded  at  the  GERESS  Site  1988-ljl89 
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The  yields  of  column  Nuttli  were  calculated  using  his  1986  paper  and  GERESS  magnitudes. 
Yields  in  columns  Ringdal  are  from  Ringdal  et  al.  (1992),  and  Ekstnxn  from  Ekstrom  and 
Richards  (1994), 
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behaviour  either  if  the  anisotropy  is  locally  very  strong  or  if  preferred  orien¬ 
tation  is  spatially  coherent  over  a  depth  range  of  several  kilometers  within 
the  topmost  part  of  the  upper  crust.  The  latter  may  for  example  be  the  case 
for  cracks  caused  by  the  regional  stress  field.  Note  that  our  P-wave  polariza¬ 
tion  data  are  sensitive  only  to  the  depth  range  where  brittle  deformation  is 
expected. 

General  anisotropic  model/hezagonal  case 

Since  the  cause  of  anisotropy  is  not  known  beforehand,  it  is  preferable 
to  find  elastic  constants  tmder  no  restricting  assumptions.  In  the  follow¬ 
ing  we  will  study  the  type  of  constraints  P-wave  polarization  data  pose  on 
elastic  constants.  We  will  restrict  our  attention  to  the  case  of  hexagonal 
symetry,  since  this  is  a  viable  model  for  most  of  the  possible  causes  of  ef¬ 
fective  anisotropy  mentioned  above.  An  exception  may  be  more  complex 
intrinsic  anisotropy  of  individual  crystals,  which  is  not  considered  here.  In 
the  case  of  a  hexagonal  or  transversely  isotropic  S3rmetry,  the  elastic  constant 
matrix  Cijkt  has  the  simple  shape 


^ijkl  — 


Cll 

Cii  —  2C66 

Cl3 

Cll  -  2C66 

Cji 

Cl3 

Cl3 

Cl3 

C33 

\  Css/ 


(7) 


with  5  unknowns,  assuming  that  the  symetry  axis  is  in  13-diiection.  Given 
this  particular  symetry,  velocity  v  and  polarization  direction  of  P-waves  de¬ 
pend  only  on  the  angle  x  between  propagation  direction  and  symetry  axis 

*3- 

3.4.1  Modelling  of  P- Velocity  and  Polarization  Patterns 
White  (1983)  and  de  Parscau  (1991b)  pve  exact  expressions  for  velocity 


pv^  =  i(cn  +  C44)sm*x  +  (^33  +  C44)cos*x+ 


(  (cil  -  C44)s»Vx  -  (<=33  -  C44)cOa*x)*  +4(c13  +  C44)*C05'X«”*X 

and  polarization  angle  (  measured  from  the  symetry  axis 

(cis  +  C4k)cosxsinx 


1/2 


tan{  = 


pv^x)  -  (cii««*X  +  C44COS*x) 


(8) 


(9) 


(8)  and  (9)  depend  only  on  the  four  elastic  constants  cn,  C13,  C33  and  C44. 
The  last  constant  c««  in  (7)  can  be  determined  only  frrom  quasi- trsmsversely 
polarized  shear  waves,  de  Parscau  (1991b)  showed  that  (9)  can  be  simplified 


to 


fan2{  = 


Vsin2x 

1  -  (1  -b  T)sin^x 


(10) 
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Figure  3.13:  Result  of  the  general  nonlinear  search:  a)  misfit  (0°to 
90°incidence),  b)  optimum  model  in  the  partial  lower  hemisphere  dis¬ 
play  (0®to  55°).  With  the  symetry  axis  at  203.3°and  41.05°.  the  (high- 
velocity)  symmetry  plane  has  sttike  of  about  113°  and  dips  with  about 
58°  to  the  N'orth.  The  optimum  set  of  parameters  is  r)  =  0.887  and 
T  =  1.09. 
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Appendix  1*1:  Chronological  Status  List  GERESS 

Dec.  1993: 

-  outage  HUB  (TelekcHii  modem  damaged  by  thundentorm) 

-  power  outage  HUB  (power  line  damaged  during  storm) 

-  power  outage  on  UPS  university  (fuse  blown  by  short  elsewhere) 

-  2  min  data  gap  possibly  due  to  power  outage  and  defective  USV  (HUB) 

-  Dec.  13  -  14,  Dec.  27  -  Jan  1.  maintenance  visits  at  HUB. 

-  data  link  GERESS  -  Kjeller  and  GERESS  -  Bochum  switched  from  land-lines 
(German  Telekom)  to  satellite  (Spaceline  Inc.). 

Jan.  1994: 

-  2  min  data  gap  possibly  due  to  power  outage  and  defective  USV  (HUB) 

-  2  thunderstorms  at  HUB  (3  times  2  min  outages  array,  station  outages) 

Feb.  1994: 

-  tenqxirary  malfunction  of  STS'2  broad-band  instrument 
Mar.  1994: 

-  2  min  data  g^s  due  to  power  outage  and  defective  USV  (HUB) 

-  implementation  of  new  amplitude-distance  curve  for  magnitude  calculation 

-  modification  of  system  kernel  on  data  acquisition  wmkstation 

Apr  1994: 

-  2  min  data  g^s  due  to  power  outage  and  defective  USV  (HUB) 

-  defective  preamplifier  C2  (broad-band  element) 

-  S  thunderstorms  at  HUB  (station  outages) 

-  Apr.  10-14  maintenance  visit  at  HUB. 

-  implementation  of  German  3.  level  seismotectonic  names  of  lASFEl  regionaliza¬ 
tion  (Leydecker  1994) 

May.  1994: 

-  3  thunderstorms  at  HUB  (array  and  station  outages) 

-  power  outage  HUB  (maintenance  by  electricity  coaq>any) 

-  May  2-6  maintenance  visit  at  HUB  (preaiiq>lifier  C2  exchanged,  calibration 
problems  fixed  on  B3,  C7,  D6,  D7,  D8,  D9,  rephuement  of  defective  USV) 

-  qxmtaneous  desynchroniziaion  on  D7 

-  temporary  malfimction  of  STS'2  broad-bond  instrument 

-  new  USV  defective  (manual  intervention  required  after  surge  or  extended  power 
outage) 
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Eksaiw  and  Richards  (1994)  estimated  source  parameters.  Considering  the 
variability  in  published  yield  determinations,  our  single  station  value  can  give  a 
valid  first  ordv  estimate  on  yield.  Hence,  for  the  2S-eleinent  GERESS  array,  a 
detection  threshold  in  the  order  of  2  kt  is  estimated. 


2.8  CONCLUSIONS 

During  the  GERESS  site  survey  (1988-1989)  and  4  1/2  years  of  GERESS 
operation  (Jan.  1990  -  June  1994),  44  presumed  nuclear  tests  have  been 
recorded:  18  from  Nevada,  14  from  the  Tuamotu  Archipelago,  6  fiom  eastern 
Kazakhstan,  5  fiom  Lop  Nor,  and  1  fiom  Novaya  Zemlya.  The  detection  thres¬ 
holds  at  GERESS  are  estimated  as  m|,  =  4.3  for  Nevada,  nO),  =  4.0  for  Tuamotu 
Archipelago,  m{,  =  4.6  fin:  Lq>  Not,  iU),  =  4.3  for  Novaya  Zemlya,  and  m{,  = 
4.5  for  eastern  Kazakhstan.  These  pafh  specific  detecticni  thresholds  deviate 
substantially  fimn  those  derived  fiom  globd  amplitude-distance  relations  (i.c., 
Gutenberg-Richter).  They  can  also  not  be  explained  with  a  simple  statiOTi  bias 
because  GERESS  generally  shows  a  negative  station  residual  of  0.6  magnitude 
units  +/-  0.35  (Harjes  et  al.,  1994).  Based  oti  2  known  yields,  a  magnitude 
yield  relation  was  calculated  for  Nevada,  indicating  that  events  with  a  yield 
estimate  in  the  order  of  1  kt  can  be  detected  at  GERESS.  Based  on  yield  esti¬ 
mates  by  DSIR  (Smith,  pers.  COTnm.),  a  magnitude  yield  relatiOTi  was  deter¬ 
mined  for  events  in  the  Tuamotu  Archipelago.  The  detection  threshold  of  mi,  = 
4.0  would  correspond  to  events  with  a  yield  estimate  in  the  order  of  0.7  kt 
However,  for  Novaya  Zemlya  and  the  Kazakhstan  test  sites,  yields  have  been 
determined  by  various  authors.  Aldiough  the  data  quantity  and  quality  at 
GERESS  for  these  test  sites  is  not  ideal,  we  estimated  a  detectioo  threshold  of 
3  kt  for  Novaya  Zemlya  and  2  kt  fOT  Kazakhstan.  Based  oti  a  cross-cofrelation 
method  and  a  cluster  analysis,  GERESS  records  can  be  used  to  discriminate 
between  the  Tuamotu  and  Fangatau&  events. 

This  case  study  also  shows  diat  caution  should  be  exercised  by  using 
results  of  theoretical  studies  on  network  detection  oqtabilities.  Instead,  path 
specific  station  detection  thresholds  should  be  used  for  a  fixed  station  netwoik 
to  derive  reliable  network  detection  ciqtabilities. 

A  reliable  predictitm  of  die  performance  of  a  future  global  monitoring  net¬ 
work  would  greatly  boiefit  fimn  information  about  exact  parameters  of  previ¬ 
ous  nuclear  tests,  fit  this  respect,  the  Chinese  and  French  authorities  should  be 
appealed  to  release  the  correspmiding  information  even  without  exact  yield 
infoimation. 


Figure  2-1;  Location  of  the  GERESS  array. 
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using  the  two  parameters 


Cl3  +  C44 

(11) 

C33  —  C44 

Cxl  —  C44 

r  - - 

(12) 

C33 — C44 

This  shows  that  we  will  i.g.  not  be  able  to  independently  determine  the 
four  elastic  constants  cn,  cia,  C33  and  C44  from  P-polarization  data.  Instead, 
constraints  on  these  constants  must  be  derived  &om  the  rj  and  r  using  the 
nonlinear  relations  (11)  and  (12).  In  addition  to  rj  and  r,  the  orientation 
of  the  symetry  axis  is  tinknown.  Therefore  we  have  2  additional  unknowns 
representing  strike  and  dip  of  the  symetry  axis.  Both  r  and  17  can  be  seen 
as  measures  of  anisotropy.  Since  cn  and  C33  give  pv^  parallel  and  normal  to 
the  symetry  planes,  r  gives  a  measure  of  the  velocity  deviation  between  these 
special  directions,  rj  is  more  difficult  to  interpret,  but  it  can  be  approximately 
related  to  the  P-velocity  anisotropy  at  small  incidence  angles  nearly  normal 
to  the  symetry  planes  (de  Parscau,  ISQlb).  For  an  isotropic  mediiim  we  have 
ri  =  T  =  l. 

A  number  of  criteria  on  the  elastic  coefficients  in  (7)  can  be  derived  from 
requiring  stability,  namely  the  condition  that  deformation  is  associated  with 
nonnegative  internal  energy.  The  requirement  is  that  (7)  is  positive  semidef- 
inite  (Backus,  1962).  It  follows  that  all  diagonal  elements  are  nonnegative. 
Furthermore,  cu  >  cesi  ^33(011  —  cce)  >  cjs*  In  our  context,  bounds  on  elas¬ 
tic  coefficients  are  important  for  finding  bounds  on  the  parameters  rj  and  r 
and  for  rejecting  physically  impossible  models.  Postma’s  (1955)  inequality 

(Cu  -  C44)(C33  —  C44)  >  (ci3  +  C44)*  (13) 

shows  that  r  >  0.  Since  r  gives  a  measure  of  the  size  of  P-wave  anisotropy,  an 
upper  bound  may  be  given  from  an  estimate  of  maximum  possible  anisotropy. 
The  sign  of  C13  +  C44,  the  nominator  of  r)  is  foimd  to  have  substantial  impor¬ 
tance  for  the  polarization  behaviour.  Cases  with  cis  +  C44  <  0  are  physically 
{>ossible,  but  they  are  fotmd  only  is  very  carefully  plmned  laboratory  ex¬ 
periments  (Helbig  and  Schoenberg,  1988).  In  these  cases,  the  polarization 
angle  is  found  to  rotate  in  a  sense  opposite  to  the  propagation  direction,  in 
contrast  to  the  normal  case  with  rotation  in  the  same  sense.  For  real  Earth 
studies  it  is  reasonable  to  assume  C13  -t-  C44  >  0,  which  gives  us  a  lower  bound 
f;  >  0.  An  upper  bound  on  17  can  be  foimd  by  applying  a  requirement  that 
C11/C44  &nd  C33/C44  be  within  a  certain  range. 

Now  we  will  illustrate  polarization  behaviour  for  two  special  cases,  namely 
(tf,  r)=  (1.43, 2.13)  and  (0.82, 1.00).  These  cases  correspond  to  elastic  param¬ 
eters  for  models  HCDl  and  HCSl  (Crampin,  1984),  which  give  the  effective 
anisotropy  associated  with  a  model  of  oriented  cracks  (Hudson,  1981)  in  the 
limiting  cases  of  fluid-saturated  and  dry  conditions.  To  derive  this  model 
of  oriented  cracks,  Hudson  (1980)  used  mean  fleld  theory  (Keller,  1964)  and 
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b)  GERESS  Array  Configuration 

622000  623000  624000  625000  626000  627000 


Figure  3.1.  a)  Geographic  location  of  the  GERESS  array  in  Southeastern 
Germany,  b)  Array  locations  in  Gauss-Kniger  coordinates.  Data  from 
25  vertical-component  shortperiod  instruments  and  three-component 
instruments  (triangles  and  rectangles)  are  used  in  this  paper.  The 
topographic  variation  across  the  array  is  about  200  meters. 
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KEA-JVE  600ni  #12 


Figure  2-9:  Map  of  event  locations  at  Nevada  test  site  (using  DOE  informa¬ 
tions,  NEIC). 
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computed  from  the  elements  of  the  covariance  matrix  axe  also  indicative  of 
the  phase  character.  The  values  of  H/V=l  and  R/T=l  are  shown  by  lines 
indicating  polarization  angles  near  45**.  For  the  S-phase  train  H/V  shows 
very  clear  horizontally  polarized  phases  5n,  Sg  and  Lg,  which  also  show 
small  values  of  R/T.  This  radial-transverse  ratio  for  P  phases  shows  clearly 
radially  polsuized  energy.  In  summary,  polarization  is  an  important  tool  for 
phase  identification.  In  the  following  we  want  to  use  polarization  data  for 
the  set  of  P  phases.  Figure  3.2  shows  that  polarization  attributes  for  P- 
phases  are  relatively  slowly  varying.  In  fact,  they  are  also  insensitive  to  the 
particular  choice  of  window  length.  On  the  6  three-component  stations  we 
observe  for  P„  a  backazimuth  of  8a  =  240.4“  ±  5.5“  and  an  incidence  angle 
of  =  44.9“  ±  4.2“ .  These  angle  estimates  correspond  closely  to  values 
extracted  manually  from  hodograms  and  are  near  predicted  values  from  the 
I.4SP91  Earth  model  and  source  locations  from  table  3.1  {8p  =  235.17“  and 
=  40.02“). 

Typically  we  have  data  from  six  three-component  stations  available.  We 
computed  standard  errors  simply  from  the  empirical  distribution  of  those 
values.  The  average  uncertainty  for  all  phases  in  the  study  is  3.5“  for  the 
incidence  angle  and  4.9“  for  the  backazimuth,  where  we  used  only  events  with 
high  Signal-to-noise  ratio. 

A  slightly  different  approach  was  given  by  Jurkewicz  (1988),  who  proposed 
averaging  covariance  matrices  of  individual  three-component  stations  before 
performing  the  eigenvector  analysis  to  stabilize  the  estimates. 

Polarization  analysis  lends  itself  to  very  efBcient  semi-automatic  process¬ 
ing.  Caution  is  required,  however,  for  phases  closely  spaced  in  time,  particu¬ 
larly  if  their  amplitudes  differ  substantially.  For  regional  phases  this  case  is 
not  imcommon.  Window  lengths  in  this  study  are  allowed  to  vary  depending 
on  data  type  and  quality.  The  median  window  length  is  near  3.5  seconds. 
Similarly,  comer  frequencies  of  a  zero-phase  bandpass  filter  were  allowed  to 
vary  within  0.6  to  4  Hz  to  adapt  to  the  different  wave  types. 

3.2.4  The  Mislocation  Hemisphere 

Figure  3.3  shows  part  of  the  lower  hemisphere  (0“  to  50“ )  of  polarization 
anomalies,  namely  deviations  of  observed  from  predicted  polarization,  where 
predicted  again  refers  to  values  computed  from  IASP91  model  with  source 
locations  from  table  3.1.  Here  the  free-surface  effect  is  taken  into  {u;coimt, 
as  will  be  discussed  below.  In  such  a  display,  the  center  corresponds  to 
vertical  incidence.  The  perimeter  gives  50“  incidence  with  North  at  the 
top.  We  choose  this  type  of  diagram  to  facilitate  comparison  with  predicted 
polarization  patterns  frxnn  anisotropic  models. 


Multiple  nudMr  events  from  two  source  repons  (Nevada,  Tuamotu)  give 
rather  OQiisuAMIt|||dAv>*tM>ts  of  polarization  to  Northern  directions.  Also 
other  ceomiiiiiflkMI  wHh  several  nearby  events  give  rise  to  consistent 

•tMfOy  incident  phases  from  the  North  there  is 
rjlorthem  directions,  while  the  Southern  part 
,  apparently  varying  rapidly  over 
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Figure  3.9:  Lower  hemisphere  representation  of  P-wave  velocity  and  polariza¬ 
tion  anomalies  for  (»?,  r)  =  (1.43,  2.13)  corresponding  to  the  cry- crack 
model  HCDl  of  Crampin  (1984). 
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Misfit  depending  on  Eta  and  Tau 
Theta=203  JOO  Phi=41.0500  eta=0^000  uu=I.09<>00  S=38^02 
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Figure  3.14:  Misfit  around  the  optimum  model  (black  square)  keeping  the 
orientation  fixed.  The  lowest  contour  line  shows  the  confidence  region 
around  the  best  fitting  model.  Nvunbers  0  to  13  refer  to  percentages  of 
mica  content  in  a  gneiss  model  (see  text). 


Shear-Wave  Velocities  for  the  Optimum  Model 
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Figure  3.16:  Prediction  of  .shear  wave  velocitie.s  for  the  optimum  model  on 
the  full  lower  hemisphere  (O'’to  90®).  qSP  is  polari/ed  i)nrall(d  to  the 
symetry  plane  and  «iSR  perpendicular. 
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